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IN  SITU  CHARACTERIZATION  OF  SOILS 
FOR  PREDICTION  OF  STRESS  STRAIN  RELATIONSHIP  OF  SOFT  CLAY 

1.  Introduction. 

The  difficulty  of  characterizing  particulate  systems  in  their  composite  form  has 
resulted  in  the  use  of  destructive  tests  to  index  the  grain  and  aggregate  properties  of 
soils  in  situ  and  to  the  arbitrary  division  of  soils  into  groups  based  on  size  and  to  the 
classification  of  soils  based  on  grain  size  distribution  and  index  tests  such  as  Atterberg 
Limits.  Due  to  the  arbitrary  nature  of  the  current  methods  of  classifying  soils,  any 
correlations  established  between  the  conventional  indices  used  to  classify  soils  and 
mechanical  behavior  need  to  be  periodically  changed  to  account  for  any  deviations. 
This  problem  can  be  overcome  if  the  method  of  classifying  soils  can  index  the 
compositional  and  aggregate  properties  of  soils  in  their  composite  form.  This  can  only 
be  achieved  by  non-destructive  means.  One  possible  way  of  achieving  the  objectives 
of 

1)  classifying  soils  in  their  composite  form 

2)  predicting  the  mechanical  properties  and  engineering  behavior  of  soils 
is  by  the  use  of  electrical  methods. 

The  objectives  of  this  report  are: 

1)  characterization  of  soil  by  electrical  methods 

2)  quantification  of  the  compositional  and  aggregate  properties  of  soils  in  terms 
of  electrical  parameters 

3)  quantification  of  the  inter  and  intra  cluster  void  ratios 

4)  establishment  of  correlations  between  electrical  parameters  and  compression 
index  X,  swell  index  K,  and  the  slope  of  the  critical  state  line,  M 

5)  methodology  to  predict  stress  strain  relationship  of  clays  in  situ. 

2,  Characterization  of  Soils  By  Electrical  Methods. 

When  an  alternating  electric  field  is  applied  to  a  clay-water-electrolyte  system, 
a  response  is  produced  which  can  be  measured  in  terms  of  a  resistance,  R,  and  a 
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capacitance,  C.  The  measured  value  of  the  capacitance  can  be  converted  into  a 


quantity  known  as  the  dielectric  constant.  This  value  is  defined  as  in  which 

Co  is  the  capacitance  of  a  condenser  with  only  a  vacuum  between  the  electrodes. 
The  dielectric  constant  is  actually  a  measure  of  the  ability  of  the  clay  to  store 
electrical  potential  energy  under  the  influence  of  an  electric  field.  From  a  knowledge 
of  the  dimensions  of  the  sample  the  dielectric  constant  and  conductivity  can  be 


calculated  from  the  following  relationships. 


e'€ 

a  A 


^-k 


(1) 

(2) 


in  which  d  =  length  of  a  specimen;  A  =  cross-sectional  area;  and  e  =  the  dielectric 
constant  of  vacuum  (8.85  x  10*^^  F/cm). 

The  dielectric  constant  of  a  dry  silicate  material  is  4.5,  and  that  of  water  is 
about  79.  A  mixture  of  soil  and  water  should,  therefore,  have  a  dielectric  constant 
weighted  between  4.5  and  79.  However,  when  the  dieJectic  constant  of  a  clay-water- 
electrolyte  system  is  measured  with  an  alternating  current  in  the  radio-frequency  range, 
it  is  found  to  be  in  excess  of  the  sum  of  the  weighted  dielectric  constant  of  the 
components.  This  measured  value,  e',  referred  to  as  the  "apparent  dielectric  constant" 
reflects  the  heterogeneous  nature  of  the  path  of  the  current  and  the  electrical  properties 
of  the  pore  fluid  and  the  clay  mineral  as  shown  by  Arulanandan  and  Smith  (1973). 

When  the  conductivity,  a,  and  apparent  dielectric  constant,  e',  of  a  cohesive 
system  are  measured  as  a  function  of  frequency,  in  the  radio  frequency  range,  it  is 
found  that  e'  and  a  vary  (see  Fig.  2).  This  variation  of  o  and  e'  is  referred  to  as 
electrical  dispersion.  However,  for  a  granular  system  it  has  been  observed  that  o  and 
e'  are  independent  of  frequency  (see  Fig.  1).  It  is  thus  possible  to  classify  all  particulate 
systems  into  two  groups.  All  cohesive  systems  characterized  by  electrical  dispersion 
behavior  and  all  granular  systems  by  nondispersive  electrical  behavior. 
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This  method  of  classification  is  both  quantitative  and  fundamental  in  nature  as 
opposed  to  the  present  qualitative  arbitrary  methods  of  classifying  soils  based  on  size 
only. 

The  factors  that  influence  the  magnitude  of  dielectric  dispersion  (see  Fig.  1) 
were  investigated  in  detail  by  Arulanandan  et  al.  (1975).  It  was  shown  that  Ae^  depend 
predominately  on  the  type  and  amount  of  clay  mineral  and  the  fabric  having  a  second 
order  effect.  Thus  the  composition  of  cohesive  soils  can  be  characterized 
nondestructively  and  quantified  in  terms  of  the  parameters  Ae^.  The  significance  of 
this  parameter  to  the  classification  of  soils  Arulanandan  et  al.  (1975),  to  the  prediction 
of  cation  exchange  capacity  Fernando  et  al.  (1975),  swell  potential  Basu  and  Arulanandan 
(1973),  lime  reactivity  Castle  and  Arulanandan  (1979),  erosion  potential  Heinzen  and 
Arulanandan  (1977)  and  compression  index  Scharlin  (1974)  has  been  demonstrated. 

3.  Electrical  Models  to  Predict  Electrical  Dispersive  Behavior 

Arulanandan  and  Mitra  (1970),  Arulanandan  and  Smith  (1973)  have  shown  it  is 
possible  to  describe  the  electrical  dispersive  behavior  of  cohesive  systems  by  the  use 
of  a  three  element  electrical  network  model. 

The  current  through  the  sample  is  considered  to  have  three  paths  (Fig.  3):  (1) 

through  the  inter  cluster  solution  and  clusters,  (2)  through  clusters  in  contact  with 
each  other,  and  (3)  through  the  inter  cluster  solution.  The  corresponding  three  element 
model  representation  of  the  impedances  of  individual  components  is  shown  in  Figure 

4. 

Referring  to  Figure  4,  k^  =  conductivity  of  the  solution,  e^.  =  dielectric  constant 
of  the  cluster,  e^  =  dielectric  constant  of  the  solution,  a  =  area  of  cluster-solution 
path,  b  =  area  of  the  particle  contact  path  and  c  =  area  of  the  solution  path.  In  the 
case  of  soUs  saturated  in  water.  Eg  «  approximately  79  and  is  approximately  equal 
to  the  conductivity  of  the  pore  fluid. 

The  apparent  dielectric  constant,  e',  and  the  apparent  conductivity  o,  for  the 
model  can  be  evaluated  by  elementary  network  analysis.  The  equations  relating  e'  and 
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a  to  the  frequency  depend  on  the  compositional  and  geometric  parameters  of  the 


system,  as  follows: 

a 


^'th  ■  d(l-d)s 


°th  d( l-d)s 


f  1  *  1 
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where 


s  = 


Hr 

1-d  d 


+■  0)  a 


2^2  I  ^  Ir 
1-d  d 


(3a) 


where 

a  =  the  capacitance  of  a  unit  capacitor  in  vacuum  (8.85  x  10”^f/cm) 
and 

u)  =  angular  frequency. 

A  digital  computer  optimization  program  utilizing  the  simplex  method  was 
employed  to  analyze  the  dielectric  constant  and  conductivity  dispersion  curves.  The 
optimization  program  determines  the  values  of  the  geometrical  (a,  b,  c,  d)  and 
compositional  (e^.,  k^,  k^,  e^)  parameters  of  the  model  shown  in  Figure  4  such  that  the 
frequency  dispersion  curves  of  the  apparent  dielectric  constant  and  conductivity 
calculated  from  the  theoretical  Equations  3,  4  fit  the  experimental  results  within 
acceptable  error  limits,  e^,  the  dielectric  constant  of  water,  was  kept  constant  at  79 
during  all  computer  runs.  The  conductivity  of  the  pore  fluid,  k^,  is  easily  determined, 
and  for  the  initial  choice  of  the  parameters,  k^  was  made  equal  to  k^.  The  initial 
value  of  "b"  is  estimated  to  be  in  the  order  of  0.01  as  "b",  the  area  of  contact,  is 
very  nearly  equal  to  zero.  The  value  of  "a"  is  dependent  on  c  since  a  +  b  +  c  =  1. 
A  typical  example  of  the  comparison  of  experimental  and  theoretical  dispersion  curves 
is  shown  in  Figure  5.  There  are,  in  fact,  six  random  parameters  c,  d,  e^,  k^,  and 
kj  involved  in  the  curve  fitting  procedure  by  optimization  technique.  At  first  sight, 
it  might  seem  that  the  choice  of  six  seemingly  random  parameters  would  be  enough 
to  fit  almost  any  set  of  curves  and  that  the  physical  meaning  of  the  result  is  doubtful. 
It  is  importrnt  to  realize,  however,  that  only  the  three  independent  geometrical 
parameters  are  genuine  unknowns  (b,  c,  d)  and  as  b  =  0,  there  are  only  two  independent 


unknown  parameters,  "c"  and  "d".  The  respective  conductivities,  and  k^,  of  the 
cluster  and  the  interstitial  solution  and  the  dielectric  constant  of  the  solid  particles 
must  have  values  that  are  reasonable  from  the  point  of  view  of  physical  science.  One 
would  expect  the  conductivity  of  the  interstitial  solution  to  be  of  the  same  order  as 
the  conductivity  of  the  pore  fluid  extract;  the  conductivity  of  the  solid  should  be  of 
the  same  order  as  the  isoconductivity  value  (Arulanandan  et  al.,  1973);  the  dieletric 
constant  of  the  clusters  should  be  between  4.5  (the  dielectric  constant  of  dry  silicate 
mineral)  and  79  (the  dielectric  constant  of  pure  water),  depending  on  the  water  content. 
Although  these  limitations  were  not  written  into  the  computer  program,  the  numerous 
results  obtained  in  earlier  investigations  show  that  the  requirements  were  indeed  met. 
However,  a  large  number  of  computer  runs  are  necessary  to  select  the  best  set  of 
parameters  that  fit  the  experimental  curves  and  are  reasonable  from  a  physical  science 
point  of  view. 

Arulanandan  (1980)  proposed  an  approximate  expression  to  quantify  the  intra 
cluster  void  ratio  ej,  utilizing  the  dielectric  constant  of  the  cluster,  e^.,  determined 
by  the  above  procedure.  Assuming  a  simple  mixing  rule,  i.e.,  that  dielectric  constant 
of  the  solution  plus  the  dielectric  constant  of  the  solids  in  a  mixture  is  equal  to  the 
dielectric  constant  of  the  mixture,  e^,  then 

Cr  =  n,  .  (1  -  n,)  (5) 

where  nj  is  the  intra  cluster  porosity.  Assuming  =  79,  ~  ^>5  (dielectric 

constant  of  dry  silicate  mineral).  Equation  5  becomes 


=  74.5  nj  +  4.5  (5a) 

Thus  by  determining  of  a  clay  water  electrolyte  system  by  use  of  the  three-element 
model  proposed  above  the  value  of  nj  can  be  quantified,  and  thus,  the  intra  cluster 
void  ratio 
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Thus  the  three-eJement  electrical  mcxlel  can  be  used  to  quantify  the  inter  and  intra 
cluster  void  ratios. 


The  quantification  of  the  inter  and  intra  cluster  void  ratios  by  the  above  approach 
is  time-consuming,  as  several  computer  runs  are  needed  to  be  made  to  obtain  the 
model  parameters,  combined  with  a  certain  amount  of  judgement  needed  to  choose  the 
right  parameters.  This  is  due  to  the  fact  that  Equations  3  and  k  are  highly  non-linear. 
In  order  to  make  the  optimization  technique  more  efficient  and  to  choose  the  suitable 
model  parameters,  it  was  thought  to  be  necessary  to  impose  an  additional  constraint 
in  the  optimization.  This  would  be  achieved  utilizing  the  measured  wate.  o.itent  of 
the  soil  as  is  explained  below.  The  intra  cluster  void  ratio,  Oj,  is  expres;  <  i  terms 
of  the  dielectric  constant  of  the  clusters,  c^,  as  given  by  Equations  5,  (  inter 

cluster  void  ratio,  Op,  is  expressed  in  terms  of  model  parameters  by  .  .^ting  the 
physical  dimensions  of  the  soil  sample  and  the  electrical  dimensions  associated  with 
the  three  current  paths  and  using  the  assumptions  made  by  Olsen  (1961).  The  total 
void  ratio  can  then  be  expressed  in  terms  of  the  three-element  model  parameters. 
The  additional  constraint  incorporated  in  the  optimization  procedure  is  as  follows: 
Measured  water  content  =  water  content  calculated  from  the  three-element  model 
parameters. 

This  constraint  overcomes  one  more  difficulty  that  could  be  encountered  as 
follows:  It  was  found  that  the  dielectric  constant  and  conductivity  curves  were  lowered 
as  the  water  content  of  the  soil  is  decreased  if  the  conductivity  of  the  pore  fluid 
remains  the  same.  However,  if  the  conductivity  of  the  soil  with  lower  water  content 
is  higher,  it  is  possible  to  obtain  identical  set  of  curves  at  both  higher  and  lower 
water  contents.  Since  the  conductivity  of  the  pore  fluid  is  allowed  to  vary  in  the 
optimization,  one  would  obtain  same  set  of  parameter  for  the  above  cases  although 
the  water  contents  are  different. 
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3a.  Inter  Cluster  Void  Ratio  as  a  Function  of  Three-Element  Model  Parameters 

Olsen  (1961)  has  used  a  three-element  electrical  model,  similar  to  the  one 
described  in  the  preceding  sections,  in  his  study  of  hydraulic  flow  through  saturated 
clays.  The  total  current  through  the  soil  sample  was  considered  to  consist  of  three 
components  as  shown  in  Figure  6.  The  derivation  of  the  expressions  for  dielectric 
constant,  e',  and  the  conductivity,  a,  are  given  in  Appendix  I.  Referring  to  Figure  6, 
=  area  of  the  solid-solution  current  path,  ^2^  =  area  of  the  contact  area  current 
psth,  A^g  =  area  of  the  solution  current  path,  =  length  of  the  current  paths  through 
each  one  of  the  elements.  The  assumptions  made  are: 

1)  The  volume  of  the  soil  sample  (A  x  L)  is  equal  to  the  volume  of  the 
solid-solution  path  (Aj^  x  L^). 

2)  The  volume  of  the  solution  path  of  solid-solution  element  (Aj^  x  ^2^ 

equal  to  the  volume  of  the  solution  path  of  the  solution  element  (a^^  x  L^). 

2  2 

Utilizing  these  assumptions  and  defining  that  a  =  1/t  ,  b  =  A2g/A|gt  ,  c  = 
A  /A  t^ 

3e  le  ,  d  =  1-d  =  '^2e^^e’  ^  expressions  for  dielectric 

constant  and  conductivity  became  same  as  the  ones  derived  by  Arulanandan  et  al. 
(1973).  However,  as  the  consequence  of  imposing  the  assumptions  made  by  Olsen,  it 

was  found  that  c  -  (l-d)a,  whereas  in  the  previous  model  a  +  b  +  c  =  1.  The  model 
parameters  are  interpreted  as  follows: 


A 

c  _  /■  .X  3e  _  Area  of  the  solution  path 
a  ■  ”  ”  Aj^  '  Area  of  the  so  lid- solution  path 

A 

b  _  2e  _  Area  of  the  contact  area  path 
a  ■  Aj^  "  Area  of  the  so  lid- solution  path 


3) 


_ Electrical  length 

Actual  length  of  sample 


=  tortuosity 


4)  Cj.  is  the  dielectric  constant  of  the  clusters 
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It  should  be  noted  that  the  tortuosity,  defined  as  t  =  L^/L,  is  an  average  value. 
This  is  due  to  the  assumption  that  the  length  of  solid-solution  current  path  is  equal 
to  the  length  of  solution  current  path.  But  one  would  expect  the  solution  current 
path  to  be  much  more  tortuous  than  the  solid  current  path. 

The  inter  cluster  void  ratio,  Op,  can  now  be  expressed  in  terms  of  the  three- 
element  electrical  model  parameters  using  Olsen's  cluster  model,  shown  in  Figure  7. 
With  the  notations  used  in  Figure  7  and  assuming  that  all  the  clusters  are  identical, 


V,/N  V,  V,  , 


I  ■  V^/N  ■  ■  V  V,  .  V  ■  1  .  e, 

C  C  s  I  s  1 


e,  = 


1  1  -  n, 


V  V 

_B _ E 


1 


p  ■  V-,  ■  V  V  V,  V  ■  1  +  e  +  e, 

p+C  s_p  l  s  p  I 

V.  V  *  V, 
s  s  s 

Q  -  — P -  ( 1  +  e  ) 

®p  -  I  -  np  ^ 


Referring  to  Figure  6 


^le'-2e  '-2e 

"p '  '  lt  '  ‘  " 
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Combining  Equations  9  and  10 

e  =  ^  (1  -  d)  (1  +  e.) 
P  d  I 


From  Equations  5a  and  6 
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‘ "  s  -  s 
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The  water  content  of  the  soil  can  be  ca  nlated  knowing  the  specific  gravity  of  the 
soil,  G^,  as  follows; 

water  contact  =  e/G^  (1^) 

The  validity  of  the  model,  when  Olsen's  assumptions  are  incorporated,  was  examined 
by  optimizing  some  experimental  electrical  dispersion  data  obtained  from  Smith,  S. 
(1971).  It  was  found  that  there  are  more  than  one  feasible  solutions  as  given  in  Table 
1.  However,  it  was  noted  that  one  of  these  feasible  solutions  predicts  a  water  content 
that  is  close  to  the  experimentally  determined  water  content  as  given  in  Table  2. 

This  observation  was  confirmed  by  optimizing  several  other  experimental 
dispersion  data  as  presented  in  Tables  3  and  4.  This  shows  that  the  model  can  predict 
the  water  content  of  the  soil. 

4.  Experimental  Procedure 

A  detailed  experimental  study  was  conducted  in  order  to  investigate  the  variation 
of  electrical  parameters  such  as  Fy,  Fj^,  A  and  the  electrical  dispersion  behavior  of 
soils  as  a  function  of  water  content.  For  this  purpose,  it  was  necessary  to  develop 
a  cell  so  that  the  above  properties  can  be  measured  as  the  soil  is  being  consolidated. 

Two  identical  square  cells  of  5"  x  1  3/8"  x  1  3/8"  (Figure  8),  fabricated  out  of 
plexiglass  of  1/2"  thickness,  were  used  to  measure  the  electrical  properties  of  the 
soils. 

One  would  be  used  for  the  determination  of  the  electrical  properties  in  the 
vertical  direction  using  two  porous  brass  plates  as  electrodes  placed  at  top  and  bottom 
of  the  sample.  The  second  cell  would  be  used  for  the  determination  of  the  electrical 
properties  in  the  horizontal  direction  using  two  brass  plates  as  electrodes  embedded 
into  the  plexiglass  walls  on  either  end  of  the  sample.  The  electrodes  in  the  cells 
would  be  connected  to  the  impedance  analyzer  by  a  6-inch  long  coaxial  transmission 
cable. 
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Since  the  impedances  of  the  transmission  line  and  the  cell  would  become  very 
significant  in  the  radio  frequency  range,  they  need  to  be  determined  accurately  and 
accounted  for  in  the  determination  of  the  dielectric  constant  and  the  conductivity  of 
the  sample.  This  is  achieved  by  representing  the  system  by  an  appropriate  equivalent 
electrical  network.  The  calibration  procedures  to  determine  the  imjjedance 
characteristics  of  the  transmission  line  and  the  cell,  using  water  which  has  a  dielectric 
constant  of  79,  is  given  in  Appendix  2.  Water  does  not  exhibit  dispersion  behavior, 
and  its  conductivity  can  be  determined  independently.  These  known  electrical  properties 
of  water  are  utilized  for  the  calibration  of  the  cells. 

A  Hewlett-Packard  impedance  analyzer  (Figure  9),  which  is  designed  for  the 
precise  impedance  measurements  in  the  frequency  range  of  1  MHz  -  1000  MHz,  was 
used  for  the  measurements  in  the  radio  frequency  range.  A  low  frequency  Impedance 
Comparator  (Figure  10)  was  used  for  the  determination  of  the  formation  factor  at  a 
frequency  of  1  KHz. 

The  soils  used  were: 

1.  Marysville  red  soil 

2.  Yolo  loam 

3.  Na-Illite 

4.  Snow  Cal  (95%)  +  Bentonite  (5%) 

5.  Snow  Cal  (60%)  +  Illite  (40%) 


The  electrical  dispersion  curves  at  different  water  contents  are  presented  in 
Figures  11-30.  It  could  be  observed  from  these  figures  thast  the  variation  of  the 
magnitude  of  dielectric  dispersion  with  water  content,  Ac^,  is  negligible.  Conductivity 
changes  considerably,  but  the  change  in  dielectric  constant  is  not  very  significant. 

The  optimized  model  parameters  using  the  experimental  dispersion  results 
corresponding  to  different  water  contents  are  presented  in  Tables  9-12.  The  anisotropy 


of  the  soil  is  mainly  reflected  in  the  model  parameter  "a"  which  represents  the 
tortuosity  of  the  current  path.  The  variations  of  the  ratios  of  the  intra  and  inter 


void  ratios,  Cj/Op,  are  almost  independent  of  the  direction  of  measurement.  This 


observation  further  confirms  the  validity  of  using  the  assumptions  made  by  Olsen  (1%1), 
since  the  void  ratios  are  scalar  quantities,  and  consequently  the  ratio,  Cj/ep,  should 
be  independent  of  the  direction  of  measurement.  The  relationship  between  ej/ep  and 
the  total  void  ratio,  e,  for  different  soils  is  shown  in  Figure  32. 

There  are  distinct  linear  relationships  between  Cj/ep  Anandarajha  (1982)  and  e 
intersecting  the  e-axis  at  a  void  ratio  of  about  0.47.  The  ratio  Cj/Op  can  only  be 
zero  when  ej  =  0.  At  this  consolidation  state,  it  would  appear  that  the  total  void 
ratio  corresponds  to  a  void  ratio  associated  with  the  maximum  packing  of  primary 
particles.  This  value  of  e  =  0.47  is  very  nearly  equal  to  the  void  ratio  associated 
with  the  maximum  packing  of  spheres  (e  =  0.43).  This  observation  provides  further 
justification  of  the  procedures  used  to  quantify  the  inter  and  intra  cluster  void  ratios. 

The  average  formation  factor,  F,  was  plotted  against  the  total  porosity,  n, 
(Tables  5-8)  as  shown  in  Figure  33.  These  results  indicate  that  there  exists  a  unique 
relationship  between  F  and  n  for  a  given  soil.  This  unique  relationships  shows  that 
the  approximate  relationship  for  sands,  F  =  n"^,  developed  using  Bruggemen’s  integration 
technique  (Dafalias  et  al.,  1979)  holds  also  for  clays.  The  results  also  show  that  the 
shape  factors,  f,  are  dependent  on  the  type  of  soils.  It  could  further  be  observed 
from  the  results  presented  in  Tables  5-8  that  the  anisotropy  index.  A,  increases  with 
decreasing  void  ratio  and  reaches  a  constant  value  at  low  stress  levels. 

5.  Electrical  Properties  in  Relation  to  Mechanical  Properties. 

It  has  been  shown  that  the  principal  factors  influencing  the  dielectric  dispersion 

6  8 

of  find  grained  soils  in  the  radio  frequency  range  (10  -  10  HZ)  are  the  compositional 
properties  of  the  different  phases  and  the  heterogeneous  nature  of  the  system.  The 
engineering  properties  of  fine  grained  soils  are  also  controlled  by  these  factors.  This 
mutual  dependency  is  used  as  a  basis  to  correlate  the  electrical  properties  to  mechanical 
properties. 
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The  purpose  of  this  research  is  to  establish  possible  correlations  between  those 
parameters  which  are  needed  for  the  prediction  of  stress-strain  behavior  of  find  grained 
soils  and  the  appropriate  electrical  parameters-  It  will  then  be  possible  to  determine 
the  stress-strain  behavior  of  fine  grained  soils  by  estimating  these  mechanical  properties 
from  the  above  correlation  and  using  them  in  a  suitable  constitutive  model.  This 
approach  facilitates  the  determination  of  the  stress-strain  behavior  of  soils  by  non¬ 
destructive  in  situ  measurements  since  the  electrical  properties  can  be  determined  in 
situ  without  disturbing  the  structure  of  the  soil. 

Any  constitutive  model  for  the  prediction  of  stress-strain  behavior  of  fine  grained 
soils  would  require  the  knowledge  of  the  slope  of  the  consolidation  line,  X,  the  slope 
of  the  swelling  line,  k,  the  slope  of  the  critical  state  line,  M,  the  initial  void  ratio, 
e^,  and  the  overconsolidation  ratio  (along  with  other  model  parameters  appropriate  for 
the  model  under  consideration).  These  parameters  are  defined  in  Figure  34.  The 
water  content  of  the  fine  grained  soils  can  be  determined  by  conventional  procedures 
and  the  initial  void  ratio,  e^,  could  be  evaluated.  The  correlations  between  the 
electrical  properties  and  the  mechanical  properties  X,  ic  and  M  for  normally  consolidated 
soils  would  be  established  in  the  following  sections. 

3a.  Compression  Index,  X 

Bolt  (1956)  has  been  able  to  predict  the  compression  characteristics  of  clays 
using  Gouy-Chapman  diffuse  double  layer  theory  and  Van't  Hoff's  theories  of  parallel 
platy  particles.  Compression  characteristics  of  Na-Montmorillonite  and  Na-illite  were 
found  to  be  close  to  the  predicted  one.  However,  the  theory  was  found  to  be  useful 
only  in  the  case  of  clays  exhibiting  very  strong  colloidal  properties  such  as 
montmorillonite.  Deviations  from  Bolt's  findings  have  been  reported  by  Mitchell  (1960) 
and  Olsen  and  Mitronovas  (1962)  and  are  ascribed  mainly  to  particle  orientation.  Quigley 
and  Thompson  (1966)  have  observed  fabric  changes  in  natural  Leda  clay  during 
consolidation  using  x-ray  diffraction  methods.  It  has  further  been  shown  by  Rosenquist 
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(1959)  that  the  compressibility  of  clay  is  dependent  on  the  type  as  well  as  the  valence 
and  concentration  of  ions  adsorbed  on  the  surface  of  the  clay  particle.  From  this 
discussion  it  is  evident  that  the  compression  index  will  depend  on  the  type  and  amount 
of  clay  mineral  and  the  electrical  interaction  between  clay  particles  and  the  resulting 
structure  of  the  soils. 

The  factors  that  influence  the  magnitude  of  dielectric  constant,  were 

investigated  in  detail  by  Arulanandan  et  al.  (1973).  It  was  shown  that  Ae^  predominantly 
depends  on  the  type  and  amount  of  clay  minera?  and  the  structure  of  soils.  The  water 
content  and  the  pore  fluid  composition  have  negligible  effect  on  the  magnitude  of 
dielectric  dispersion.  The  experimental  results  presented  in  Figures  11-30  further 
confirm  this  finding.  Therefore  one  would  expect  a  correlation  between  the  slope  of 
isotropic  consolidation  line,  X  and  the  magnitude  of  dielectric  dispersion,  because 

of  their  mutual  dependency.  Sharlin,  J.R.  (1972)  has  shown  that  such  relationship 
indeed  exists  in  the  case  of  fine  grained  soils  as  shown  in  Figure  35.  The  relationship 
between  X,  obtained  from  the  isotropic  consolidation  tests  for  the  soils  reported  in 
this  research  (Figure  36)  and  Ae^  of  the  respective  soils  shows  a  close  agreement  with 
the  earlier  finding  of  Sharlin,  3.F.  (1971).  This  confirms  the  general  validity  of  the 
relationship  between  the  slope  of  isotropic  line,  X  and  the  magnitude  of  dielectric 
dispersion,  Ae^ 

5b.  Swell  Index,  k 

Lambe  (1960)  discussed  the  physico-chemical  forces  acting  on  a  clay  particle. 
The  fluid  pressure  between  two  clay  particles  is  higher  than  the  pore  pressure  of  the 
free  water  due  to  the  relatively  higher  cation  concentration  between  the  particles. 
The  pressure  difference  between  two  points,  one  within  two  clay  particles  and  the 
other  in  the  free  water,  is  considered  to  be  numerically  equal  to  the  electrical  repulsive 
forces  acting  between  two  platy  clay  particles.  If  these  repulsive  forces  are  very 
strong  in  a  given  type  of  soil,  one  would  expect  this  soil  to  swell  considerably  on  the 
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removal  of  the  external  loads  so  that  the  equilibrium  forces  acting  on  a  particle  would 
be  maintained.  The  voids  between  two  clay  particles  would  be  larger,  if  these  repulsive 
forces  are  strong,  compared  to  the  inter  particle  voids  in  soils  with  weaker  repulsive 
forces. 


Thus,  with  the  cluster  concept  (Michael  et  al.,  1954  and  Quirk,  1959)  one  could 
argue  that  the  soil  which  has  higher  proportions  of  intra  cluster  void  ratio  with  respect 
to  the  inter  cluster  void  ratio  at  a  particular  total  void  ratio,  would  swell  more  on 
removal  of  load  compared  to  the  soil  with  lower  proportions.  The  ratio  of  intra  and 
inter  cluster  void  ratios,  Cj/Op  can  not  be  quantified  considering  two  parallel  particles 
since  the  clay  particles  aggregate  and  exist  in  a  random  arrangement.  The  electrical 
method  of  quantifying  ®|/Cp  accounts  for  the  structure  of  the  soils,  since  no  assumptions 
are  made  regarding  the  arrangement  of  clay  particles  in  quantifying  the  inter  and  intra 
cluster  void  ratios  using  the  electrical  methods. 

On  the  basis  of  the  preceding  discussion,  the  swell  index,  k,  was  correlated 
with  Cj/ep  at  a  common  total  void  ratio  of  e  =  1.0,  i.e.  (ej/ep)^  _  j  q.  The  values 
of  e,/e  at  e  =  1.0  for  different  soils  were  obtained  from  Figure  32  and  the  values 

1  p 

of  K  were  obtained  from  Figure  36.  A  linear  relationship  between  Cj/Op  at  e  =  1.0 
and  K  is  evident  on  a  log-log  space  as  shown  in  Figure  37  Anandarajha  (1982). 


5c.  Slope  of  the  Critical  State  Line,  M 

When  a  soil  element  is  sheared  under  drained  or  undrained  conditions,  experimental 
results  indicate  that  the  soil  element  fails  when  the  stress  path  reaches  the  critical 
state  line  independent  of  the  initial  stress  state  of  the  soil  element  (Schofield  et  al., 
1968)  as  illustrated  in  Figure  34.  At  failure,  the  void  ratio  and  the  effective  mean 
normal  pressure  lie  on  a  unique  line,  referred  to  as  a  critical  state  line,  on  the  e-p 
space  as  shown  in  Figure  34.  This  concept,  widely  known  as  critical  state  concept 
(Roscoe  et  al.,  1968)  is  the  basis  of  the  Cam-clay  theory  developed  by  the  Cambridge 
group.  From  this  discussion,  it  is  evident  that  the  slope  of  the  critical  state  line,  on 
the  p-q  space,  represents  the  ultimate  shear  strength  of  the  soils. 
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Lambe  (1960)  has  discussed  the  factors  controlling  the  shear  strength  of  soils. 
The  major  component  of  ultimate  shear  resistance  of  fine  grained  soil  is  considered 
to  be  friction.  Figure  37  illustrates  different  ways  in  which  dilatancy  and  frictional 
behavior  could  be  mobilized.  One  could  therefore  say  that  "M"  is  a  function  of  particle 
size,  shape,  surface  texture  and  the  structure  of  soils  resulting  from  the  electrical 
attractive  and  repulsive  forces  between  clay  particles. 

The  shape  factor,  T,  defined  in  the  expression,  F  =  n'^,  (Dafalias  et  al.,  1979) 
is  a  function  of  the  shape  of  the  particles  and  cementation  in  the  case  of  sands. 
However,  in  clays,  the  shape  factor  would  reflect  the  physico-chemical  interaction 
between  particles  which  depends  on  shape,  composition  and  arrangement  of  particles. 

The  factors  influencing  the  formation  factor  can  also  be  examined  from  the 
theoretical  equations  for  conductivity  obtained  from  the  three-element  electrical  model. 
From  Equations  3a  and  4, 


0)  =  0  -  kjil  -  Op)  +  np 


+  bk^  +  ck^ 


(1  -  Hp)  .  (j^)  np 


bnp  >  b(j^)2  np  >  anp  (1  -  *  anp  2(^)  ^  a(^) 


(15) 


(15a) 


2  2 
where  a  =  1/t  ,  np  =  1  -  d  =  inter  cluster  porosity,  bt^  =  cluster  contact  area/area 

of  the  solution-cluster  path. 

It  is  approximately  true  that  F  =  _  q-  Therefore, 

k  k 

F  =  - -  =  function  (n  ,  t,  b)  (16) 

The  ratio  of  the  conductivity  of  the  cluster  to  the  conductivity  of  solution  k^/k^, 
reflects  the  intra  cluster  void  ratio  and  fabric  of  the  soil  and,  t,  the  tortuosity  reflects 
the  pcirticle  shape,  size  and  the  fabric  of  the  soil. 


-16- 


Therefore,  it  could  be  argued  that  the  shape  factor,  f,  which  relates  the  formation 
factor  to  the  total  porosity  is  a  function  of  the  shape,  size  and  texture  of  the  clay 
particles,  and  the  fabric  of  the  soils.  The  emisotropy  index.  A,  was  found  to  be  another 
parameter  controlling  the  strength  in  the  case  of  sands.  Soils  which  exhibit  higher 
anisotropy  would  have  lower  strength  when  measurements  are  made  in  the  direction 
of  the  major  principal  axis. 

Based  on  these  considerations,  it  is  considered  that  an  electrical  index,  A  /T, 
may  be  related  to  "M",  since  the  factors  influencing  these  two  parameters  are  the 
same. 

In  order  to  be  able  to  predict  the  stress  path  in  the  p-q  space  and  the  stress-strain 
behavior  in  the  q-e  space,  it  is  necessary  to  obtain  "M"  so  that  when  values  of  X  k, 
and  M  are  used  in  a  chosen  constitutive  model,  the  predicted  and  the  observed 
stress-strain  relationships  would  be  reasonably  close.  Two  constitutive  models,  namely 
the  basic  Cam-clay  model  developed  by  the  Cambridge  group  and  a  bounding  surface 
plasticity  model  developed  by  Dafalias  et  al.  (1980)  were  chosen  for  this  purpose. 

The  procedures  recommended  by  Schofield  et  al.  (1968)  were  followed  in  obtaining 
the  "M"  values.  According  to  Cam-clay  theory,  the  relationship  between  Alog^p  and 
q/p  should  be  linear,  where  the  stress  variables  correspond  to  undrained  conditions. 
An  example  of  the  experimental  relationship  between  Xlog^p  and  q/p  for  the  soil  type 
Snow  Cal  (60%)  +  lllite  (40%)  is  shown  in  Figure  39.  The  slope  of  this  straight  line 
is  M(X  -  ic).  Thus,  knowing  X  and  k,  the  slope  of  critical  state  line,  M,  can  be 
determined. 

A  calibration  procedure  (Herrmann  et  a!.,  1981)  was  followed  in  obtaining  "M" 
values  using  bounding  surface  plasticity  model.  The  details  of  these  theories  would 
be  found  in  the  latter  part  of  this  report. 

2 

The  relationship  between  the  electrical  index,  A  /T,  (Table  13)  and  the  slope  of 
the  critical  state  line,  M,  obtained  from  Cam-clay  theory  and  the  bounding  surface 
plasticity  theory  are  shown  in  Figures  40  and  41,  respectively.  The  reasonably  good 
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2  _ 

non-linear  relationship  beteen  M  and  A  /f  indicates  that  the  factors  controlling  "M", 
which  is  a  mechanical  property,  and  A  /T,  which  is  an  electrical  property,  are  related 
to  each  other. 

6.  In  Situ  Prediction  of  Stress-Strain  Behavior  of  Cohesive  Soils 

Application  of  the  theory  of  plasticity  has  gained  wide  popularity  over  the  past 
few  years  for  the  analytical  prediction  of  stress-strain  relationship  of  soils.  The  most 
notable  and  earliest  of  the  constitutive  models  developed  using  the  theory  of  plasticity 
is  the  Cam-clay  theory  developed  by  the  Cambridge  group  (Schofield  et  al.,  1968). 
The  Cam-clay  model  can  represent  the  strain  hardening  behavior  depending  on  the 
state  of  the  soil  with  respect  to  the  critical  state  line.  The  verification  of  the 
Cam-clay  model  under  simple  shear  and  triaxial  stress  states  has  been  carried  out 
satisfactorily  on  soils  that  are  "wet"  of  critical  state  (Roscoe  et  al.,  1968).  However, 
this  model  has  not  been  demonstrated  to  be  representative  of  the  strain  softening 
behavior  of  the  soils.  Yield  surface,  along  with  the  normality  rule,  gives  rise  to 
prediction  of  unrealistic  behaviors  off  isotropic  loading.  This  is  due  to  the  fact  that 
the  yield  surface  derived  as  the  result  of  assumptions  made  in  Cam-clay  theory  does 
not  intersect  the  p-axis  at  right  angle. 

Roscoe  et  al.  (1968)  modified  the  theory  by  introducing  an  ellipse  for  the  yield 
surface  and  using  normality  condition.  The  size  of  this  elliptic  yield  surface  is 
completely  defined  by  the  initial  isotropic  consolidation  pressure,  p^,  and  the  slope  of 
the  critical  state  line,  M.  The  value  of  R,  which  is  defined  as  the  ratio  of  the  mean 
normal  pressure,  p^,  at  which  the  yield  surface  intersects  the  p-axis  and  the  mean 
normal  pressure,  Pj,  at  which  the  yield  surface  intersects  the  critical  state  line  (refer 
to  Figure  42),  is  assumed  to  be  2.0,  a  constant  value.  It  is,  however,  desirable  to  be 
able  to  use  different  values  of  R  in  the  theory  in  order  to  be  able  to  predict  the 
stress-strain  behavior  of  a  broad  range  of  soil  types. 
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Numerous  theories  have  since  been  developed  to  describe  the  stress-strain 
behavior  of  soils  under  a  general  loading  condition.  Bounding  surface  plasticity  model 
(Dafalias  et  al.,  1980),  cap  model  (Di  Maggio  et  al.,  1971),  endochronic  theory  (Valanis 
et  al.,  1971),  and  model  developed  by  Prevost  (1978)  are  some  of  them. 

The  bounding  surface  plasticity  model  has  been  developed  to  describe  a 
generalized  stress-strain  behavior  of  fine  grained  soils  under  a  variety  of  conditions 
such  as  normally  consolidated  and  over  consolidated,  strain  hardening  and  strain 
softening,  monotonic  and  cyclic,  compression  and  extension,  loading  paths  in  the  tension 
zone  (i.e.,  the  effective  mean  normal  stress  is  tension),  etc.  The  verification  of  this 
theory  under  some  of  the  above  conditions  has  been  carried  out  satisfactorily.  The 
value  of  R  can  be  varied  as  required  to  predict  a  given  experimental  stress-strain 
relationship,  unlike  in  the  case  of  critical  state  theory  developed  by  Roscoe  et  al. 
(1968).  When  R  =  2,  these  two  theories  predict  identical  stress-strain  behavior  in  the 
case  of  normally  consolidated  soils. 

For  the  purpose  of  this  research,  the  feasibility  of  using  the  basic  Cam-clay 
theory  (Schofield  et  al.,  1968)  and  a  bounding  surface  plasticity  theory  developed  by 
Dafalias  et  al.  (1980)  for  the  prediction  of  the  stress-strain  behavior  of  normally 
consolidated  fine  grained  soils,  will  be  discussed  in  the  following  sections. 

6a.  Basic  Cam-Clay  Theory 

The  Cam-clay  theory  was  developed  utilizing  the  critical  state  concept,  as 
discussed  earlier.  The  stress  and  strain  variables  in  the  triaxial  space  are  defined  as 
follows: 
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where 


Oj  :i  total  vertical  stress 

0^  =  total  lateral  stress 

u  =  pore  pressure 

Oj'  =  effective  vertical  stress 

o^'  =  effecitve  lateral  stress 

q  =  deviatoric  stress  variable 

p  =  effective  mean  normal  pressure 

• 

e  =  deviatoric  strain  increment 

v/v  =  volumetric  strain  increment 

=  elastic  deviatoric  strain  increment  and 

=  plastic  deviatoric  strain  increment 

p  =  volumetric  stress  increment 

q*  =  deviatoric  stress  increment 

The  stress-strain  theory  has  been  derived  using  a  power  law  and  a  stability 
criteria  and  assuming  that  the  soil  does  not  display  any  recoverable  shear  strain,  as 
follows: 

Energy  balance  criteria, 

p(^)P  +  q(e)P  =  M  p  le*^!  (18) 

Normality  concept, 

p  +  q(e)^  =0  (19) 

Using  Equations  18  and  19,  the  equation  o.  the  yield  surface  can  be  derived  as  follows: 

3  =  M  logg(^)  (20) 
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where  =  precondolidation  pressure  as  defined  in  Figure  42.  The  material  behavior 
is  elastic  if  the  stress  state  lies  within  the  yield  surface  and  is  plastic  if  the  stress 
state  lies  outside.  The  relationships  between  q,  p  and  v  for  undrained  condition  are 
shown  in  Figure  42. 

The  expression  representing  the  normally  consolidated  undrained  stress-path  can 
be  derived  based  on  the  fact  that  the  volume  of  the  sample  does  not  change  during 
the  test.  This  is  given  by  Equation  (21) 


a  . 

P 


(21) 


The  undrained  stress-strain  relationship  is  given  by  Equation  (22) 


where 


•  1  X  -  K  £  T]_ 

M-n  1+e  p'^MP 


n  =  q/p 
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Using  Equat ion  21 , 


q^  =  M  Po  exp 


(22) 

(22a) 

(22b) 


(23) 


where  q^  is  the  deviatoric  stress  of  normally  consolidated  soil  at  failure. 

The  undrained  stress-strain  behavior  of  a  normally  consolidated  soil  under  triaxial 
loading  conditions  is  fully  defined  by  Equations  21  and  22.  The  information  required 
for  the  analytical  prediction  of  stress-strain  behavior  is  X,  ic ,  M,  e^  and  p^. 
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6b.  The  BoundiriR  Surface  Plasticity  Model  for  Co'.iesive  Soils 

In  this  section  the  basic  features  of  the  bounding  surface  plasticity  model  for 
cohesive  soils  will  be  presented,  emphasizing  the  concepts  rather  than  the  analytical 
formulation  which  can  be  found  in  given  references.  Particular  attention  will  be  given 
in  describing  the  role  of  the  constitutive  constants  since  these  are  the  final  objectives 
of  the  effort  to  correlate  mechanical  with  electrical  properties. 

6b.  1.  The  General  Concept 

Let  us  first  present  the  salient  features  of  the  bounding  surface  concept,  and 
in  order  to  be  specific  the  presentation  will  be  associated  with  the  particular  formulation 
for  cohesive  soils.  Assuming  plastic  isotropy,  the  analytical  expression  for  the  bounding 
surface  in  the  space  of  stress  invariants  is  given  by 

F(T,  3,  a,  e")  =  0  (24) 

where  I,  J  and  a  stand  for  the  first  stress  invariant,  the  second  deviatoric  stress 
invariant  and  the  Lode  angle  (depends  on  the  ratio  of  the  second  and  third  deviatoric 
stress  invariants),  respectively,  a  bar  over  stress  quantities  imply  stress  states  on  F=0 
and  finally  e"  represents  the  plastic  change  of  the  void  ratio  (the  change  of  the  void 
ratio  due  to  plastic  volumetric  strains).  Exact  definitions  of  the  above  can  be  found 
in  [Dafalias  et  al.  1982].  The  e"  is  the  only  plastic  internal  variable  which  controls 
the  hardening/softening  response.  A  typical  schematic  representation  of  the  bounding 
surface  is  shown  in  Fig.  43  by  means  of  2  ellipses  and  one  hyperbola  in  the  meridional 
section  of  the  I,  3  space,  but  notice  that  a  enters  also  the  description  assuming  a 
fixed  value  for  each  such  meridional  section,  on  which  the  model  parameters  N,  A 
and  R  depend  (to  be  explained  subsequently).  The  influence  of  e"  appears  through  the 
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dependence  of  1^  on  it,  being  the  intersection  with  the  positive  I  axis.  Such  a 
surface  could  also  play  the  role  of  a  classical  yield  surface  and  here  is  where  the 
concept  of  the  bounding  surface  is  introduced:  instead  of  viewing  the  surface  described 
by  eq.  (24)  as  a  locus  of  points  in  stress  space  which  sharply  delineate  tlie  purely 
elastic  (inside  the  surface)  from  the  elasto-plastic  domain  (on  the  surface),  it  is  rather 
viewed  as  a  bounding  envelope  of  elasto-plastic  stress  states  in  or  on  the  surface,  with 
an  additional  property:  the  closer  is  an  actual  stress  state  to  the  surface,  the  more 
extensive  is  the  corresponding  increment  of  plastic  deformation  for  a  given  stress 
increment. 

In  order  to  present  a  well  defined  analytical  form  of  the  above  properties,  the 
following  concepts  have  been  developed  [Dafalias  1981,  Dafalias  et  al.  1982].  For 
each  actual  stress  state  I,  3,  a  in  or  on  the  surface  (observe  that  for  a  general 

development  not  restricted  to  isotropy,  the  stress  state  should  be  defined  by  means 

of  the  stress  components  not  just  their  invariants),  Fig.  43,  a  corresponding  unique 
"image"  stress  state  1,  3,  a  on  F  =  0  is  defined  by  means  of  a  properly  defined  mapping 
rule.  In  the  present  formulation  the  "radial  mapping"  rule  is  employed  which  is  shown 
in  Fig.  43  and  which  essentially  defines  the  "image"  stress  state  by  projecting  radially 
the  current  state  I,  3,  ot  onto  F  =  0  using  a  properly  defined  stress  point  as  the  center 
of  projection.  Here,  such  a  point  will  be  chosen  to  lie  always  on  the  I-axis  and 
given  by  =  CI^  (0  _<  C  _<  1).  In  Fig.  43  the  value  C  =  0  is  used,  i.e.  the  projection 

center  is  the  stress  origin  0  itself.  Observe  that  on  the  basis  of  the  above  definition 

of  the  mapping  rule  a  =  a  since  they  both  refer  to  the  same  meridional  section.  For 
a  discussion  of  the  general  conditions  that  any  mapping  rule  must  satisfy,  the  reader 
is  referred  to  [Dafalias,  1981].  It  suffices  here  to  say  that  when  the  actual  stress 
point  lies  on  F  =  0,  its  "image"  is  identical  to  it.  Obviously  this  is  satisfied  by  the 
present  "radial  mapping"  rule. 
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The  definition  of  the  "image"  stress  serves  two  purposes.  First,  at  the  "image" 
stress  the  corresponding  gradient  \7F  defines  the  loading-unloading  direction  at  the 
actual  stress.  In  order  to  visualize  this,  one  may  think  of  the  "vector"  VF  translated 
parallel  to  itself  from  f,  3,  a  or  F  =  0  where  it  is  defined,  to  1,  3,  a  where  it  acts, 
Fig.  43.  Second,  the  Euclidean  distance  6  between  T,  3,  a  and  I,  3,  a  is  a  measure 
of  the  proximity  of  the  actual  stress  to  the  bounding  surface  and  its  effects  on  the 
plastic  stiffness  is  expressed  by; 


K  =  K  +  h(a)  H{1,  3,  a,  e") 


s6> 


(25) 


where:  Kp  =  The  actual  plastic  modulus. 

Kp  =  A  bounding  plastic  modulus  associated  with  the  "image"  stress  and 
obtained  from  the  consistency  condition  F  =  0. 
h  =  a  shape  hardening  parameter;  a  function  of  a. 

H  -a  function  of  the  state  (used  to  normalize  the  plastic  modulus), 

r  =  the  Euclidean  distance  between  the  projection  center  and  the 
"image"  stress. 

<  >  =  the  McCalley  brackets  defining  the  operation  <A>  =  AifA  > 
0  and  <A>  =  0ifA_£0. 

s  =  the  elastic  factor,  such  that  for  6  >  r/s  <r  s6  >  =  0  *  K  =  «. 

P 

Indeed  the  definition  of  F  =  0,  eq.  (24),  and  the  corresponding  definition  of  the  mapping 
rule  together  with  eq.  (25)  are  the  cardinal  elements  of  a  bounding  surface  formulation. 
For  6  =  0  (point  on  the  bounding  surface),  eq.  (25)  yields  Kp  =  Kp,  i.e.  the  bounding 
surface  behaves  like  a  yield  surface.  The  effect  of  each  of  the  above  parameters, 
will  be  considered  in  the  following  subsections  where  all  material  and  model  constants 
are  discussed  systematically. 


-24- 


6b.2.  State  Variables  and  Constitutive  Constants 

Frequently  there  is  confusion  in  the  literature  of  constitutive  modeling  on  the 
meaning  attributed  to  the  words  material  or  model  constants,  parameters  and  variables. 
The  ultimate  goal  of  a  constitutive  model  for  history  dependent  inelastic  deformation 
is  to  relate  the  stress  increment  to  the  strain  increment  (or  rates)  and  vice-versa.  To 
this  end  a  number  of  other  variables  enter  the  formulation  which  effect  the  above 
stress-strain  rate  relations  in  each  step,  embodying  basically  the  effect  of  past  loading 
history.  These  are  frequently  called  internal  variables  and  can  also  be  called  model 
parameters  in  the  sense  that  they  must  be  updated  in  each  step  in  order  to  provide 
the  proper  stress-strain  rate  relation.  These  model  parameters  are  not  constants  and 
their  number  effects  greatly  the  storage  requirement  and  efficiency  of  a  constitutive 
model.  In  the  present  formulation  the  bounding  surface  model  has  only  one  internal 
variable  or  model  parameter,  the  e"  or  respectively  the  1^  (eq.  (24)),  recalling  that 
the  changes  of  I^  and  e"  are  related.  It  follows  that  the  state  of  the  material  is 
defined  by  means  of  the  stress  and  the  internal  variables,  here  the  I^.  The  initial 
values  of  the  state  variables  must  be  known  before  any  loading  process,  but  let  us 
carefully  avoid  confusing  these  initial  values  with  the  constitutive  constants. 

The  constitutive  constants,  as  their  name  implies,  are  constants  whose  values 
are  determined  during  the  calibration  process  and  remain  constant  (for  the  given 
material)  throughout  any  subsequent  stress-strain  history.  Again  their  number  and  ease 
of  calibration  affect  the  efficiency  of  the  model.  Note  that  for  one  and  the  same 
material  the  same  values  of  these  constants  must  be  used,  inhomogeneity  due  to  the 
state  (not  the  material)  can  be  described  by  means  of  the  different  values  of  the  state 
variables  (stress  state  and  internal  variables).  Of  course,  more  general  inhomogeneity 
implies  different  material  from  location  to  location,  thus  yielding  different  constitutive 
constants  in  addition  to  possible  differences  of  the  state.  Expecting,  however  that 
large  portions  of  the  structure  under  analysis  are  made  from  the  same  material,  only 
one  set  of  constitutive  constants  must  be  stored  for  each  portion,  while  as  many 
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sets  of  state  variables  or  model  parameters  must  be  stored  and  updated,  as  is  the 
number  of  points  used  in  each  portion  for  a  finite  element  analysis.  Thus,  having  a 
small  number  of  model  parameters  is  much  more  important  for  numerical  computations 
than  a  small  number  of  constitutive  constants. 

6b.3.  Constitutive  Constants  for  the  Bounding  Surface  Model 

The  constitutive  constants  for  the  bounding  surface  model  in  its  present  form 
can  be  divided  in  3  distinct  categories. 

6b. 3a.  The  classical  material  constants 

These  are  constants  related  to  a  specific  material  property  which  is  invariant 
with  respect  to  loading  history,  within  the  limits  of  approximation  and  experimental 
error.  As  such,  they  can  be  determined  by  classical  experiments. 

The  first  is  the  consolidation  index  A,  i.e.  the  slope  of  the  nomal  consolidation 
line  in  the  e-ln(p)  plot. 

The  second  is  the  rebound  or  swell  index  <,  i.e.  the  slope  of  the  swell  and 
recompression  line  in  the  e-ln(p)  plot,  and  is  directly  related  to  the  elastic  bulk  modulus 
K  for  elastic  strain  increments.  The  A  and  k  define  in  combination  the  change  of 
with  respect  to  changes  of  e",  and  more  precisely  the  above  relation  depends  on  their 
difference  A  -  This  is  expected  since  A  is  a  "tangent"  modulus  encompassing  both 

plastic  and  elastic  deformations,  while  k  is  only  reflecting  elastic  volumetric  strains, 
thus  A  -  K  naturally  appears  related  to  the  change  of  plastic  strains  only  and  therefore 
to  the  hardening  process  (change  of  with  respect  to  e"). 

The  third  constant  is  the  elastic  shear  modulus  G  which  can  be  obtained  from 
the  initial  slope  of  a  deviatoric  stress  versus  deviatoric  strain  curve  divided  by  three, 
in  a  typical  triaxial  experiment.  An  alternative  way  is  to  consider  instead  the  third 
constant  to  be  Poisson's  ratio  v  and  obtain  G  indirectly  by  mean  of  v  and  the  bulk 
modulus  K  (depending  on  p  and  k).  This  yields  a  variable  G  with  p,  which  although 
giving  better  results  in  many  cases  does  not  satisfy  the  requirement  that  elastic 
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strain  should  be  obtained  from  a  potential  of  the  stresses  and  vice-versa.  The  reason 
for  this  apparent  inconsistency  is  that  K  and  G  are  more  complex  function  of  the 
stress  state  than  presently  assumed. 

The  fourth  material  constant  is  the  critical  r‘ate  line  slope  (compression) 
in  q-p  space,  related  to  the  friction  angle.  Although  a  corresponding  (extension) 
can  be  obtained  from  by  assuming  Coulomb  criterion,  it  is  more  accurate  to  specify 
as  an  independent  material  constant  by  extension  experiments.  Observe  that  in 
the  invariant  stress  space  N  =  M/3/3,  Fig.  43. 

6b.3b.  Model  constants  associated  with  the  shape  of  .he  bounding  surface 

The  model  constant  considered  in  this  and  the  following  subsection,  influence 

the  material  response  indirectly  and  are  not  associated  with  any  distinct  scalar  material 

property,  thus  they  are  called  model  constants  rather  than  material  constants.  Fig. 

43  will  clarify  the  contents  of  this  subsection. 

The  shape  of  the  bounding  surface  must  satisfy  certain  general  requirements 

such  as  it  must  intersect  the  I  axis  at  points  (I  and  1  )  where  the  nomal  to  the  surface 

o  t 

is  along  the  axis  (isotropy),  and  at  the  intersection  with  the  C5L  (Critical  State  Line) 
the  nomal  must  be  parallel  to  the  J  axis  (critical  failure  with  zero  volumetric  strain). 
Besides  these  two  requirements,  there  is  no  definite  reason  to  prefer  one  shape  from 
another.  In  the  course  of  the  development  of  this  model  it  was  found  that  one  shape 
which  can  provide  a  very  good  fit  of  experimental  data  was  obtained  by  a  combination 
of  2  ellipses  and  one  hyperbola  as  shown  in  Fig.  43. 

Ellipse  //I  is  defined  by  the  variable  R,  function  of  a,  by  means  of  its  values 
in  compression  and  R^  in  extension  and  an  interpolation  rule  for  all  other  values 
of  a.  The  R  determines  the  Ij  =  I^/R  and  hereby  the  shape  of  the  ellipse,  since  1^-lj 
and  3j  are  the  two  principal  axes.  The  values  of  R^,  and  R^  can  be  obtained  by 
curve  fitting  the  equation  of  the  undrained  stress  path  in  triaxial  compression  and 
extension,  respectively.  Notice  also  that  R  measures  the  "distance"  between  the  nomal 
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consolidation  and  the  critical  state  lines  in  the  e-lnp  plot.  Ellipse  al  is  closely  related 
to  the  behavior  of  nomally  consolidated  or  slightly  overconsolidated  samples  in  the  wet 
side  of  the  CSL,  as  is  obvious  from  its  position  with  respect  to 

If  one  now  attempts  to  use  the  extension  of  Ellipse  al  to  the  left  of  the  point 
C,  Fig.  U3,  for  heavily  overconsolidated  samples,  it  will  be  impossible  to  properly  fit 
experimental  data  for  this  range.  Experience  shows  that  a  surface  closer  to  the  CSL 
is  required,  and  as  such  a  hyperbola  is  introduced  with  an  asymptote  parallel  to  CSL. 
Its  proximity  to  CSL  is  controlled  by  the  scalar  quantity  A,  Fig.  43,  which  in  general 
depends  on  a  and,  therefore,  requires  two  model  constants  for  its  determination,  i.e. 
its  values  A^,  A^  in  compression/extension.  These  model  constants  can  be  determined 
by  fitting  the  stress-strain  curves  at  large  OCR,  especially  by  fitting  the  peak  stress. 
The  smaller  is  the  values  of  A,  the  lower  is  the  value  of  peak  stress  etc. 

Finally,  the  bounding  surface  extends  into  the  purely  tensile  range  (I<0, 
compressive  stress:  positive)  as  Ellipse  //2,  requiring  one  more  model  constant  T  which 
defines  its  intersection  I^  with  the  1  axis  by  1^  =  TI^.  There  is  not  much  experimental 
evidence  to  be  used  for  the  accurate  determination  of  T  at  this  stage,  but  it  is 
expected  that  T  is  proportional  to  the  cohesive  strength  of  the  material. 

In  summary,  the  following  5  model  constants  define  the  shape  of  the  bounding 
surface;  R^,  R^,  A^,  A^  and  T.  It  is  expected  that  the  use  of  only  three  constants 
R  =  (R^  +  Rg)/2,  A  =  (A^  +  Ag)/2  and  T  can  give  satisfactory  agreement  with 
experimental  data.  Moreover,  often  A  and  T  may  be  fixed  at  typical  values  of  0.02 
and  0.1,  respectively,  leaving  only  R  as  the  only  important  model  constant  for  the 
shape  of  the  surface. 
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6b. 3c.  Model  constants  associated  with  response  for  OCR  >  1 

The  first  such  constant  is  C  defining  the  projection  center  1^  =  Cl^.  In  Fig. 
43  the  particular  case  for  C  =  0  is  only  shown.  The  determination  of  C  requires  at 
least  two  undrained  stress  paths  at  small  and  large  OCR,  since  the  position  of  1^ 
essencially  determines  the  shape  of  these  paths  in  the  q  -  p  (or  3  -  I)  space.  From 
experience,  it  is  expected  to  be  between  0.3  and  0.5. 

The  remaining  constants  are  associated  with  eq.  (25).  The  value  of  s>l,  defines 
a  range  of  purely  elastic  response  within  the  bounding  surface,  called  the  elastic  nucleus 
since  for  6>r/s  one  has  ®  according  to  the  definition  of  <  >.  Again  its  value 

can  be  obtained  by  observing  the  shape  of  the  undrained  stress  path,  namely  when  it 
deviates  from  a  straight-up  line  in  the  p-q  space  it  indicates  the  boundary  of  the 
elastic  nucleus.  Its  value  controls  also  cylic  stabilization.  It  is  p,  jibJe  to  fix  the 
value  of  s  at  a  typical  value  between  1.5  and  2. 

And  finally,  2  additional  and  perhaps  mot  important  constants  for  OCR  >  1  define 
the  values  of  the  shape  hardening  parameter  h  in  compression  (h^)  and  extension  (h^), 
thus  h  is  obtained  by  interpolation.  It  is  again  expected  that  a  single  constant  h  = 
(h^  +  h^)/2  will  provide  very  good  results.  The  greater  is  the  value  of  h,  the  stiffer 
is  the  response  within  F  =  0,  and  as  a  matter  of  fact  as  h  -►  «>  one  retrieves  the 
classical  yield  surface  formulation  since  then  from  eq,  (25)  for  points  within 

the  bounding  surface. 

In  summary,  the  following  4  constants  are  necessary  for  OCR  >  1  :  C,  s,  h^ 

and  h^. 
e 


6b.3d.  A  final  conclusion 

Averaging  the  values  of  different  constants  in  compression  and  extension  and 
fixing  some  of  them  at  typical  values,  one  finally  ends  up  with  the  following  set  of 
9  material  and  model  constants:  X,  k  ,  G,  M,  R,  A,  C,  s  and  h. 

Even  further.  A,  C  and  s  can  be  fixed  for  classes  of  typical  clay  soils,  and  one 
is  left  with  only  6  constants;  X,  k,  G,  M,  R  and  h.  The  effort  of  this  project  is  to 
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obtain  a  correlation  between  electrical  measurements  and  the  above  constants,  thus 
enabling  us  to  predict  the  soil  constitutive  model  by  in-situ  measurements  (the  electrical 
ones)  and  hence  be  able  to  perform  analysis  of  earth  structures  on  the  basis  of  the 
calibrated  bounding  surface  soil  plasticity  model. 

6b.4.  Laboratory  Calibration  of  the  Bounding  Surface  Model 

In  order  to  provide  checks  for  the  electrical  measurement  calibration  procedure 
that  is  being  developed  as  part  of  this  research,  means  must  be  available  for  calibrating 
the  model  using  standard  laboratory,  triaxial  test  data.  In  this  section  a  brief  description 
of  a  computer  aided  semi-automated  calibration  scheme  for  using  such  data  is  given; 
greater  detail  and  examples  are  given  in  reference  (DeNatale,  et  al.  1983,  DeNatale, 
1982).  If  different  behavior  in  extension  and  compression  is  to  be  modeled  the  scheme 
requires  that  data  from  at  least  six  triaxial  tests  (three  in  extension  and  three  in 
compression)  be  available,  in  addition  to  isotropic  consolidation  and  rebound  data.  If 
similar  (or  average)  behavior  in  compression  and  extension  is  assumed  then  only  three 
test  results  need  be  available. 

In  formulations  such  as  the  Bounding  Surface  model,  which  employ  a  small 
number  of  constitutive  constant  whose  roles  in  the  constitutive  formulation  are  each 
well  defined,  the  calibration  process  becomes  systematic  and  straightforward.  However, 
reliance  on  user  expertise  is  still  high,  since  all  manual  curve  fitting  procedures,  by 
their  very  nature,  require  both  judegement  (in  deciding  just  what  constitutes  the  "best" 
overall  fit)  and  familiarity  (in  deciding  how  much  each  constants  value  must  be  changed 
to  improve  a  given  prediction). 

In  order  to  simplify  the  model  calibration  process,  a  computer  code  has  recently 
been  developed  by  DeNatale  (1982)  and  tested  on  a  variety  of  real  soils.  The  code 
employs  a  Quasi-Newton  optimization  strategy  to  locate  that  set  of  paramter  values 
which  minimizes  the  discrepancy  between  the  model  predictions  and  the  experimental 
observations  included  in  the  calibration  data  base.  Because  this  new  computer  aided 
prcedure  greatly  reduces  the  dependence  of  calibration  success  on  user  expertise,  it 
significantly  increases  the  accessibility  and  usefulness  of  sophisticated  material  models 
to  the  general  engineering  community. 


6b.4a.  Objective  function 

Since  the  calibration  of  a  material  model  involves  minimizing  the  error,  or 
residual,  between  the  observed  and  predicted  soil  response,  the  process  can  quite 
naturally  be  viewed  as  an  optimization  problem.  Therefore,  in  order  to  develop  a 
computer  directed  calibration  procedure,  it  is  necessary  to  (i)  construct  an  objective 
finction  to  serve  as  a  scalar  measure  of  tl«  goodness  of  a  particular  solution,  and  (ii) 
select  a  search  strategy  to  enable  the  minimum  of  this  function  to  be  located  in  an 
efficient  and  reliable  manner. 

For  those  cases  in  which  all  aspects  of  the  soil  response  can  be  expressed  in 
terms  of  a  single  quantity  (by  using,  for  example,  q  vs  Cj,  u  vs  Cj  and  p'  vs  Cj  rather 
than  q  vs  p',  q  vs  and  u  vs  Cj),  a  vertical  measure  of  the  residual  may  be  reasonable. 
However,  if  the  distinction  between  independent  and  dependent  variables  cannot  be 
made  clear,  there  is  no  more  reason  to  use  a  vertical  measure  (yry(x))  than  there  is 
to  use  a  horizontal  measure  (x=x(y)).  In  these  cases  it  is  probably  more  appropriate 
to  use  a  measure  such  as  the  shortest  distance  between  the  experimental  observation 
and  prediction  curve.  While  there  is  no  reason  to  suspect  that  this  alternative  is 
theoretically  more  sound,  such  a  "Euclidean"  measure  is  probably  closer  to  what  one 
intuitively  uses  when  estimating  the  error  between  two  curves. 

The  calibration  code  permits  either  absolute  or  squared  residuals,  to  be  used  to 
form  the  objective  function.  When  these  two  options  are  combined  with  the  choice 
of  either  a  vertical  or  a  Euclidean  measure  (as  described  above),  the  user  has  four 
possible  ways  to  define  the  error  at  a  point.  Recent  research  by  DeNatale  (1982)  has 
shown  that  the  location  of  the  global  minimum  remains  essentially  the  same,  regardless 
of  which  options  are  used.  However,  preliminary  applications  to  a  variety  of  aritifical 
and  real  soils  indicated  that  the  absolute  Euclidean  measure  results  in  a  more  well- 
behaved  objective  function  that  can  be  most  easily  minimized. 
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In  order  to  combine  the  residuals  at  various  points  whithin  a  given  relation  (for 


example,  the  q  vs  relation  from  a  test  at  OCR=I),  or  from  relations  of  the  same 
kind  obtained  from  different  tests  (for  example,  the  u  vs  £j  relation  from  tests  at 
OCR=l  and  2),  it  is  necessary  to  first  define  what  is  meant  by  "equal"  error.  A  given 
solution  could  be  defined  to  be  equally  good  at  two  points  A  and  B  if  either  (i)  there 
was  the  same  absolute  error  at  both  A  and  B,  or  (ii)  there  was  the  same  relative  error 
at  A  as  at  B.  The  decision  as  to  which  definition  of  equality  should  be  used  is  entirely 
up  to  the  user.  The  computer  code  permits  the  use  of  either  extreme,  or  any  strategy 
in  between.  Consequences  of  the  various  choices  are  discussed  further  in  DeNatale 
(1982).  It  should  be  noted  that  most  physical  analogies  lie  somewhere  in  the  middle 
-  a  given  dial  gauge  or  pore  pressure  transducer  may  be  accurate  to  within  a9b,  but 
may  fail  to  register  any  meaningful  readings  below  a  magnitude  of  g. 

The  code  permits  any  number  of  tests,  relations  and/or  individual  observations 
to  be  included  in  the  calibration  data  base.  Because  the  various  response  relations 
will  generally  be  of  different  dimensions  (such  as  stress  vs  stress,  stress  vs  strain, 
strain  vs  strain,  etc.),  all  data  is  nondimensionalized  so  that  the  errors  from  different 
relations  can  be  properly  combined.  Different  weights  may  be  assigned  to  specific 
components  of  the  data  base  if  it  is  felt  that  certain  tests,  relations  or  observations 
are  more  reliable  or  representative  than  others,  or  if  it  is  necessary  to  have  the  final 
model  predictions  fit  some  data  more  closely  than  others.  The  consequences  and  proper 
role  of  weighting  factors  is  again  discussed  in  DeNatale  (1982). 

6b.4b.  Optimization  strategy 

An  extremely  large  number  of  optimization  strategies  have  been  suggested  over 
the  last  30  years,  with  the  performance  of  a  given  approach  being  strongly  dependent 
on  the  particular  type  of  problem  to  which  it  is  applied.  Hence,  it  is  generally  agreed 
that  there  is  no  single  best  algorithm,  but,  rather  only  strategies  which  perform  best 
when  applied  to  certain  classes  of  problems. 

In  selecting  the  most  suitable  approach,  a  key  factor  is  whether  or  not  first 


and  second  derivative  information  can  be  readily  obtained.  With  the  Bounding  Surface 
model,  as  with  most  sophisticated  material  models,  the  governing  equations  are  so 
complex  that  it  is  essentially  impossible  to  directly  relate  the  relevant  repsonse 
parameters  (such  as  p,  q,  u,  e^,  etc.)  to  the  constitutive  constants  employed  by  the 
formulation.  Thus,  the  objective  function  must  be  formed  by  summing  a  series  of 
discrete  weighted  residuals,  and  therefore  first  and  second  derivative  information  is 
not  explicitly  available. 

Under  these  conditions,  the  current  consensus  among  those  most  active  in 
optimization  research  is  that  a  Quasi-Newton  strategy  with  finite  difference 
approximations  to  derivatives  will,  if  properly  implemented,  generally  exhibit  the  most 
efficient  and  reliable  performance.  Hence,  a  Quasi-Newton  strategy  was  incorporated 
into  the  calibration  code  to  direct  the  search. 

Finally,  the  code  permits  constraints  to  be  imposed  on  the  various  constitutive 
constants  in  the  form  of  simple  bounds 

<  Xj  <  u, 

where  and  u.  represent  the  minimum  and  maximum  values  that  parameter  x^  can 
assume.  A  restriction  of  this  kind  would  be  appropriate  if  either  (i)  there  were  certain 
theoretical  restrictions  placed  on  the  value  which  a  given  parameter  could  assume,  or 
(ii)  there  were  certain  ranges  of  parameter  values  beyond  which  the  numerical 
implementation  of  the  material  model  became  unstable,  or  (iii)  certain  material 
properties  were  observed  experimentally  to  vary  over  some  finite  range,  and  there  was 
no  overwhelming  reason  to  fix  them  prior  to  calibration  at  any  particular  values. 

6b.  4c  Application  to  the  Bounding  Surface  Model 

The  search  for  the  optimal  set  of  constitutive  constant  is  directed  by  the 
Quasi-Newton  strategy  previously  described.  However,  in  order  to  evaluate  the  objective 
function  at  some  location  x,  it  is  necessary  to  first  generate  a  corresponding  set  of 
model  predictions.  To  accomplish  this,  the  calibration  code  relies  on  two  subroutines 
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—  EVAL  and  CLAY  —  developed  by  Herrmann  et  al.  (1980,  1982)  during  previous 
research  with  the  Bounding  Surface  model  and  briefly  described  in  the  following  section. 
Subroutine  performs,  essentially,  single  element  incrementaliterative  finite  element 
analyses  of  bodies  under  a  homogeneous  state  of  stress  and/or  strain.  Subroutine  CLAY 
consists  of  a  numerical  implementation  of  the  governing  constitutive  equations,  and 
thus,  when  called,  provides  the  appropriate  material  response  to  the  given  stress  and/or 
strain  increment. 

To  verify  the  viability  of  the  new  computer  aided  calibration  procedure,  the 
method  was  applied  to  a  number  of  representative  soil  data  bases,  both  artifical  and 
real.  The  outcome  of  these  studies  is  discussed  by  DeNatale  (1982). 

6b.5  Numerical  Implementation  of  the  Model 

As  noted  in  the  last  section,  calibrations  of  a  material  model  to  fit  laboratory 
data  requires  that  means  be  available  for  numerically  evaluating  the  model,  i.e.,  means 
for  using  the  model  to  predict  response  to  various  test  conditions.  Furthermore  if  a 
material  model  is  to  be  of  any  practical  value  numerical  means  must  be  available  for 
using  the  model  in  analyzing  engineering  structures.  To  these  ends  a  numerical  analysis 
procedure  has  been  developed  for  evaluating  the  bounding  surface  plasticity  model  for 
cohesive  soils.  This  development  is  briefly  discussed  in  the  following  sub-sections. 

6b.5a  Incremental  Stress-Strain  Relation 

The  theoretical  formulations  previously  discussed  yield  a  relationship  between 
the  rates  of  the  stress  and  strain  components  (Hermann  et  al.,  1980);  expressing  this 
relation  in  matrix  form  gives: 

{ct}  =[D]  {e}  (26) 

where  [D]  is  a  matrix  whose  dependence  on  the  stress  state,  the  internal  variable  and 
the  constitutive  constants  follow  simply  from  the  basic  equations  of  the  bounding 
surface  model  Herrmann  et  al.  (1980). 


In  general  the  response  of  an  elasto-plastic  body  is  highly  nonlinear  and  path 
dependent.  Thus  a  general  numerical  analysis  procedure  for  elasto-plasticity  problems 
requires  an  incremental  solution  and  unless  the  increments  are  made  excessively  small, 
iteration  must  be  conducted  in  each  increment  to  account  for  the  nonlinear  behavior. 

To  be  able  to  use  eq.  26  in  an  incremental  solution  procedure,  it  must  be 
expressed  in  an  incremental  form.  Consider  the  n*^  step  of  an  incremental  analysis; 
i.e.,  the  solution  has  been  found  at  n-1,  and  it  is  desired  to  calculate  the  incremental 
change  that  will  give  the  solution  at  n.  Because  of  the  nonlinear  behavior,  iteration 
is  required  to  establish  the  incremental  change.  In  the  k-1  iteration  of  this  process, 
the  estimates  of  the  stress  and  strain  states  at  n  are  given  by 


{a} 

=  lo} 

+  {Ao} 

(27) 

n  ,k-  1 

n- 1 

n,k- 1 

{t) 

=  {e} 

+  {Ae} 

(28) 

n,k-l 

n-1 

n,k-l 

Even  though  rate  independent  Isehavior  is  being  considered,  it  is  convenient  to 
think  in  terms  of  the  time  history  of  the  quantities  involved.  Integrating  eq.  (26) 
from  time  t^  ^  to  t^  and  using  the  trapezoidal  formula  to  approximate  the  right  hand 
side  gives 

{Ao}  =  [D]  {Ao}  (29) 

n,k  n,k-l  n,k 

where 

[5]  =  I  [D]  +  [D]  (30) 

n,k-l  n-1  n,k-l 

Eq.  (29)  is  the  desired  incremental  stress-strain  equation  for  iteration  k  of 
increment  n. 

Because  [Dl^  is  a  function  of  the  stress  and  strain  states  at  n,  it  is  necessary 
to  base  its  value  on  the  estimates  of  the  previous  iteration  (eqs.  (27)  and  (28)).  The 
resulting  value,  denoted  by  j  do  not  in  fact  satisfy  eq.  (30).  In  general,  the 

estimates  of  the  stress  and  strain  increments  used  in  the  calcuation  of  j 


do 


not  satisfy  eq.  (30).  This  inconsistency  disappears  as  iteration  (at  the  global  level) 
continues  and  convergence  occurs.  However,  to  speed  up  the  global  convergence 
process,  iteration  can  be  introduced  in  the  calculation  of  [D]  ,  ,  (e.g.,  for  finite 

element  applications,  iteration  when  the  properties  for  each  integration  point  are 
calculated)  to  eliminate  this  inconsistency;  in  light  of  the  additional  iteration  at  the 
system  level  a  rather  loose  convergence  criterion  appears  to  be  justified.  Thus,  the 
stress  estimate  is  iteratively  modified  (in  this  process  {Ae}  ,  ,  is  held  constant  in 
order  to  preserve  compatibility  of  the  global  strain  field)  to  give  a  refine  I  incremental 
stress  estimate  { Aa}*  .  and  a  corresponding  consistent  incremental  properties  estimate 
[D]*|^  j  which  satisfy  eq.  (29). 

The  developed  incremental  stress-strain  relationship  has  been  used  in  evaluating 
simple  laboratory  test  results  and  research  is  currently  being  carried  out  to  use  it  in 
finite  element  analysis  of  engineering  structures;  preliminary  results  indicate  that  the 
model  is  exceedingly  easy  to  implement  in  finite  element  analysis  and  is  numerically 
very  well  behaved.  The  application  to  homogeneous  tests  is  described  in  the  next 
sub-section. 

6b. 5b.  Application  to  Homogeneous  Tests 

As  previously  noted,  the  assessment  of  the  characteristics  of  a  material  model, 
and  in  the  fitting  of  it  to  experimental  measurements,  a  means  must  be  available  for 
using  the  model  to  predict  the  results  of  simple  homogeneous  tests.  For  cohesive 
soils,  program  EVAL  has  been  written  for  this  purpose.  EVAL  can  be  used  to  predict 
the  behavior  of  homogeneous  soil  samples  subjected  to  arbitrary  homogeneous  stress 
and  strain  histories  for  either  drained  or  undrained  conditions. 

The  solution  history  is  broken  into  "history  segments."  Within  each  history 

segment,  a  consistent  combination  of  six  stress  and  strain  components  are  prescribed 

(i.e.,  the  histories  of  e  or  o  ,  e  or  cr  and  y  or  t  ).  A  different  combination 

X  X  y  y  yz  yz 

may  be  prescribed  in  each  history  segment.  For  example,  a  uniaxial  test  might  involve 
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two  segments  with  a  specified  value  of  axial  strain  achieved  at  the  end  of  the  first 
segment  and  with  unloading  to  zero  axial  stress  specified  in  the  second.  Each  history 
segment  is  broken  into  increments  with  the  iteration  conducted  within  each  increment. 
The  analysis  conducted  by  EVAL  is  essentially  a  one-element  finite  element  analysis 
of  a  homogeneous  body.  The  analysis  is  straightforward  and  well  behaved  with  only 
two  special  features  worth  noting. 

Because  the  analysis  can  be  used  for  the  extreme  case  when  either  all  the  stress 
or  all  the  strain  components  are  specified,  both  the  stress  and  strain  vectors  are 
checked  for  convergence.  The  convergence  check  on  the  stress  increment,  however, 
must  be  done  with  some  care.  The  problem  is  that  the  relative  measure  of  error, 
L,  (Ao  . -Ao  ,  ,)/L,(Aa  .  ),  used  for  the  check  is  meaningless  when  applied  for 

1  n^K  rijK”!  1  rijK 

near  failure  conditions  (because  Ao.  =  0,  Lj  (Ao^l^)  =  0).  To  avoid  this  difficulty, 

the  denominator  of  the  error  measure  is  limited  to  a  minimum  value  of  L,(o  ,)/10. 

1  n- 1 

In  the  above,  the  Lj  norm  is  the  sum  of  absolute  values. 

The  second  feature  involves  the  starting  estimate  for  the  strain  increment.  For 
the  first  iteration  of  the  first  increment  of  all  history  segments  beyond  the  first, 
special  precautions  are  used  to  avoid  a  problem  of  non-uniqueness  that  may  occur  for 
unstable  soil  behavior.  With  strain-controlled  loading  it  is  possible  to  predict  soil 
response  for  the  falling  portion  of  a  stress-strain  curve.  Having  arrived  at  some  point 
in  the  unstable  region,  at  the  end  of  a  particular  history  segment,  a  user  might  switch 
to  stress-control  to  unload  the  sample  in  the  next  history  segment.  However,  starting 
from  a  point  on  a  descending  stress-strain  curve  (or  its  equivalence  for  a  multi-axial 
event)  the  strain  may  either  increase  (loading)  or  decrease  (unloading)  for  a  specified 
reduction  in  stress;  thus,  the  subsequent  behavior  is  non-unique.  Because  such  a  stress 
controlled  specification  (from  a  point  on  a  descending  stress-strain  curve)  is  only 
meaningful  for  an  unloading  situation,  the  analysis  must  select  this  path.  The  selection 
of  the  unloading  path  is  assured  by  using  for  the  initial  estimates  of  the  components 
of  the  strain  increment,  for  the  first  iteration  of  the  first  increment  of  each  new 
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history  segment,  the  negative  (reduced  by  the  multiple  .01  in  order  to  reckon  with  the 
stiffer  behavior  encountered  for  unloading)  of  the  calculated  values  for  the  last 
increment  of  the  previous  history  segment. 

As  previously  noted,  program  EVAL  is  an  integral  part  of  the  recently  developed 
model  calibration  algorithm;  it  is  used  for  evaluating  the  residuals  (differences  between 
modal  predictions  and  experimental  observations)  needed  to  define  the  objective  function 
in  the  optimization  procedure.  It  has  also  been  used  extensively  in  verification  studies 
where  model  predictions  are  compared  to  experimental  observations,  e.g.  see  Herrmann 
et  al.,  1981. 

6c.  Application  to  In  Situ  Prediction  of  Stress-Strain  Behavior 

The  values  of  the  slope  of  the  critical  state  line  were  obtained  both  from 

Cam-clay  theory  and  the  bounding  surface  plasticity  theory  and  correlated  with  an 

electrical  index  as  described  earlier.  It  can  be  observed  from  Equation  (21)  that  a 
plot  of  q/p  vs.  Xlog^p  would  be  linear  and  the  slope  of  the  line  is  M/(X  -  <).  The 
variables  q  and  p  are  the  experimentally  determined  values  during  the  undrained  test. 
Knowing  X  and  <,  the  slope  of  the  critical  state  line,  M,  can  then  be  computed.  An 

example  of  a  plot  of  Xlog^p  vs.  q/p  is  shown  in  Figure  39.  The  slope  is  2S.0  and  X 

=  0.09918  and  k  :=  0.016399  for  this  soil.  From  this  the  "M"  value  is  computed  to 
be  2.318.  The  values  of  "M"  plotted  in  Figure  40  were  obtained  in  this  manner. 

The  model  parameters  defining  the  bounding  surface  plasticity  model  can  be 
obtained  using  a  set  of  triaxial  test  data  by  a  trial  and  error  procedure,  i.e.,  the  model 
parameters  should  be  chosen  so  that  the  predicted  and  measured  stress- strain 
relationships  are  close  to  each  other.  A  computer  program  was  developed  Herrmann 
et  al.  (1982)  as  described  in  the  earlier  section  incorporating  the  bounding  surface 
theory  in  order  to  predict  the  stresses  and  strains.  This  program  was  used  to  obtain 
the  value  of  M  which  predicts  the  relationships  of  q  vs.  Ej  and  q  vs.  p  (Ej  =  axial 


strain)  as  close  to  t'-ie  corresponding  experimental  relationships  as  possible.  The  values 

2 

of  "M"  obtained  in  this  manner  are  correlated  with  the  electrical  index,  A^/T  as  shown 
in  Figure  41. 

It  could  be  observed  from  Figures  40  and  41  that  the  values  of  "M",  required 
by  these  theories  to  predict  the  experimental  stress-strain  behavior  to  best  of  their 
ability,  are  different.  This  could  be  ascribed  mainly  to  the  differences  in  the  shape 
and  size  of  the  yield  surface  in  the  Cam-clay  theory  and  the  bounding  surface  in  the 
bounding  surface  theory.  These  surfaces  are  used  as  plastic  potential  to  define  the 
direction  of  the  plastic  strain  increment  vectors  assuming  the  associated  flow  rule,  as 
discussed  below. 

From  Equations  23  and  26,  the  undrained  strength  of  the  soil  depends  on  X,  tc, 
M,  R,  and  the  initial  effective  confining  pressure,  p^.  Since  R  =  2.72  was  used  in 
the  bounding  surface  theory,  the  undrained  strengths  predicted  by  Cam-clay  theory 
would  be  different  from  the  ones  predicted  by  the  bounding  surface  theory  by  the 
ratios  of  "M".  The  undrained  strengths  normalized  with  respect  to  the  initial  effective 
confining  pressures  are  presented  in  Table  14. 

In  order  to  see  how  realistic  these  predictions  are,  the  stress-strain  relationships 
predicted  by  these  theories  were  compared  with  the  experimental  relationships 
Anandarajha  (1982)  as  shown  in  Figure  44  for  the  soil  type  Snow  Cal  (60%)  +  Illite 
(40%).  The  comparisons  are  made  on  the  q  -  Ej  space  and  the  q  -  p  space.  It  could 
be  noted  from  Figure  44  that  the  overall  deviations  of  these  predictions  from  the 
experimental  results  are  larger  in  the  case  of  Cam-clay  theory. 

From  this,  the  bounding  surface  plasticity  theor'  :'eveloped  by  Dafalias  et  al. 
(1980)  is  considered  to  be  desirable  because  of  its  realistic  predictions  and  the  generality 
of  the  theory  for  use  under  various  conditions. 

Comparison  of  the  experimental  stress-strain  relationships  and  the  bounding 


surface  predictions  for  three  other  soils,  Snow  Cal  (95%)  +  Bentonite  (5%),  Marysville 
Red  soil,  and  Illite  are  shown  in  Figures  45-47. 


A  reasonably  good  agreement  between  the  experimental  and  theoretical 
predictions  of  stress-strain  behavior  of  fine  grained  soils  is  seen  from  the  results 
presented  in  Figures  44-47.  It  is  therefore  possible  to  predict  the  stress-strain  behavior 
of  normally  consolidated  fine  grained  soils  using  bounding  surface  plasticity  model  if 
model  parameters  can  be  determined  by  in  situ  tests.  In  situ  model  parameters  could 
be  predicted  from  in  situ  electrical  measurements.  Therefore  it  can  be  concluded  that 
if  in  situ  electrical  measurements  (i.e.,  apparent  dielectric  constants  and  conductivity 
in  the  radio  frequency  range  along  with  conductivity  of  the  solution  extract)  are  made 
in  situ,  stress-strain  behavior  of  normally  consolidated  fine  grained  soils  could  be 
predicted. 

7.  Current  Work  in  Progress 

The  work  reported  above  were  for  the  prediction  of  stress  strain  behavior  of 
fine  grained  soils.  Current  study  is  concerned  with  the  stress  strain  behavior  of  mixed 
soils  (a  mixture  of  sands,  silts,  natural  and  artificial  clays  up  to  10  percent). 

It  should  also  be  noted  that  the  study  presented  in  this  report  is  limited  to  the 
prediction  of  stress  strain  behavior  of  normally  consolidated  soils.  Prediction  of  the 
in  situ  stress  strain  behavior  of  over-consolidated  soils  needs  to  be  investigated. 

8.  Future  Work 

It  is  proposed  to  investigate  the  possibility  of  the  prediction  of  initial  in  situ 

stress  state  using  the  electrical  method.  The  prediction  of  in  situ  stress-strain  behavior 

of  over  consolidated  soils  requires  a  knowledge  of  O.C.R.  or  preconsolidation  pressure 

(P^)  and  the  value  of  K.  It  is  proposed  to  establish  correlations  between  the  mechanical 

parameters  K^,  P^  and  K  and  electrical  parameters  in  order  to  enable  the  prediction 

of  K  ,  P  and  K  in  situ, 
o’  c 


L 
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9.  Summary  and  Conclusions 

The  feasibility  of  quantifying  the  inter  and  intra  cluster  void  ratios  by 
characterizing  the  soils  by  electrical  methods  is  investigated.  The  mechanical  properties 
of  fine  grained  soils  related  to  the  stress-strain  behavior  are  quantified  by  a  fundamental 
characterization  of  soils  by  electrical  methods.  The  Cam-clay  constitutive  model  and 
a  bounding  surface  model  were  used  to  predict  the  stress-strain  behavior  of  the  normally 
consolidated  fine  grained  soils.  The  bounding  surface  theory  was  found  to  yield  more 
realistic  predictions  compared  to  the  Cam-clay  theory.  This  study  shows  that: 

1.  The  inter  and  intra  cluster  void  ratios  can  be  quantified  by  characterizing  the 
soils  by  electrical  methods. 

2.  The  mechanical  properties  required  to  predict  the  stress-strain  behavior  of  soils 
can  be  determined  by  electrical  methods. 

3.  Stress-strain  behavior  of  normally  consolidated  fine  grained  soils  can  be 
determined  by  electrical  methods. 

4.  Since  the  electrical  properties  can  be  measured  in  situ,  stress-strain  behavior 
of  normally  consolidated  fine  grained  soils  can  be  determined  by  the  in  situ 
measurements  of  electrical  properties. 
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APPENDIX  1 


DERIVATION  OF  THE  PROPOSED  THREE  ELEMENT 
ELECTRICAL  MODEL  DISPERSION  EQUATIONS 

Three  current  paths  considered  and  the  associated  tfiree  elements  aie 
described  in  Fit>  2.S  in  the  text.  The  equivalent  admittance  of  each  ele-nent 
is  derived  belo\>.'. 
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Solid  path:  Resistance  R 
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Element  2 


FoJlo\^ing  the  same  procedure  as  in  element  I, 
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Now  sjnre  1  -  I|  ♦  1^  * 
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Assumptions 

1)  Assuming  that  the  volume  of  soil  sample  -  volume  of  element  1  of  elccirica 
rriodel; 
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2)  Assuming  that  the  volume  of  solution  part  of  element  1  =  volunie  of 

element  3; 
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Observe  that  the  expressions  for  c*  ,  o  ,  and  S  are  the  sa^ne  as  thos'> 

theot.’  tficor.  ^  ^  u  ... 

gi\’en  by  equations  2.30,  2.31,  and  2.32  in  the  rnain  text,  with 


a  r  b  : 


*^2e 


A 


2  ’ 


c  - 


3c 


■  le 


-  d, 


'Ze 


^  t  t 

le 


A,  t^ 
le 


l-d  and  f  -  • 


From  these. 
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A,  f 
le 
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APPENDIX  2 


CALIBRATION  OF  ELF.CTRICAL  CELLS 


Representation  of  tiie  systcni  b>  an  equivalent  elec  trical  network: 


Transmission  Line: 


L  -  distributed  inductance  per  unit  length 
R  =  distributed  resistance  per  unit  length 
C  =  stray  capacitance  per  unit  length 
G  =  stray  conductance  per  unit  length 
=  characteristic  impedance  of  tfie  line 
y  -  propagation  constant  of  the  line 
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Cell  without  Sample: 


-3  " 
'•'zi  ^ 


Cell  with  Sample: 
P 


Y5  r 


By  knowing  Zj,  Z2.  characteristic  impedance  of  the  transmission  line,  Z^^, 
propagation  constant  of  the  line,  y,  impedances  of  the  cells  and  Z^.  ti,'.- 
resistance  and  the  capacitance  of  the  sample  (Z^)  can  be  reduced  frorr.  the 
measured  value  of  the  impedance  at  a  required  frequency.  They  are  determined 
as  follou's: 

1.  The  values  of  Rp  Lj.  drid  are  given  in  the  impedance  analy?!  r 
operating  manual. 

2.  Charar tens!  1C  impedance  of  ttic  transmission  line: 

1  .  .  1 


^  "V  “open  "“short 


Y  .  i  tanh'  7 
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3. 


Impedance  of  the  cell: 


A  O' 


7 

A.  ‘ 

1 

A  O 


■  R 


“i  '  ^ 

_  tanh 

I 


Yl  -  Z 


I 


Thus  by  measuring  the  impedances  with  CC  short  and  open,  without 
sample,  and  accounting  for  the  impedance  of  the  impedance  analyzer,  cell 
impedances  2^  and  can  be  computed  using  the  above  expression. 


Cell  without  the  sample  whth  CC  short: 

1 


B 

o- 


L 


O 

B 


-or 


Cell  without  the  sample  with  CC  open: 


Now  b>  ttioasurini'  the  im|,H'd.iiir(>s  with  t'lt-  sample  in  the  cel',  resi  !  im  • 
and  the  capacitance  of  the  sample  can  be  computed  since  the  imped, trues  e; 
the  rest  of  the  system  are  known. 

Calibration  results: 

The  procedure  described  in  the  preceding  section  was  followed  iii 
determining  the  iir.pedance  characteristics  of  tlic  system.  This  procedure,  alorig 
with  the  known  values  of  the  system  impedance  parameters  w-ere  incorpo-ateo 
in  a  computer  program,  w'ritten  in  ftASIC,  wfiich  computes  the  dielectric  constant 
and  the  conduct  vity  of  the  sample  when  the  length  of  the  sample,  measur  ng 
frequency,  resistive  part  of  the  measured  impedance  and  the  reactive  part  of 
the  measured  impedance  are  used  as  input  parameters.  The  cell  impedance  was 
determined  at  different  heights  of  The  samples.  A  listing  of  the  program  is 
attached  at  the  end  of  Appendix  2. 

Water  with  known  values  of  conductivity  was  used  in  The  cells,  and  the 
measurements  were  made  in  order  to  check  the  validity  of  ttic  equivalent 
electrical  network  representation  of  the  transmission  line  and  the  cell  and  the 
accuracy  of  the  measurement  of  the  parameters  concerned.  The  water  is  found 
to  have  a  dielectric  constant  of  79,  and  it  is  not  supposed  to  exhibit  dispersion 
in  the  radio  frequency  range.  Rut  the  measured  values  of  dielectric  constant 
and  conductivities  w'ere  found  to  vary  slightly  with  the  frequency.  This  was 
suspected  to  be  due  to  the  inaccurate  determination  of  the  impedances  with  the 
electrodes  CC  short,  since  the  electrodes  had  to  be  shorted  using  a  length  of 
wire,  and  this  has  a  certain  unknown  value  of  inductance. 

Utilizing  the  fact  that  the  wat<?r  does  not  exhibit  dispersion  and  its 
electrical  properties  are  known,  cell  impedance  values  were  adjusted  slightly  bs 
trial  and  error  to  obtain  the  expected  results.  A  typical  example  of  the  computed 
dielectric  constant  and  conductivity  of  water  is  given  in  Table  A.l. 
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The  determined  citaracter istic  impedance  and  the  propagation  factor  of 
the  transmission  line  with  frequency  are  sliown  in  Figures  A.  I  and  A. 2.  The 
inductance  and  the  capacitance  of  the  cells  lor  5  different  heights  of  the  sainple 
are  presented  in  Table  A. 2.  The  resistance  of  the  cell  is  plotted  against  the 
frequency  for  different  heights  of  the  satnple  as  shown  in  Figure  A. 3.  The 
typical  example  of  the  behavior  of  the  conductance  with  frequency  is  shovi.n  in 
Figure  ,\A. 
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Table  A.l 


COMPUTED  DIELECTRIC  CONSTANT  AND  CONDUCTIVITY  OF  WATER; 
HEIGHT  OF  sample  =  2.564  inch. 

CONDUCTIVITY  OF  WATER  AT  100  KHZ  -  5.05  x  10“’  mho/cm 


F  requency 
(MHZ) 

Measured 

impedance 

Computed 
dielectric  constant 
ol  water 

Computed 
Conductivity 
of  water 

X  10'^ 

(mho/cm) 

R(R) 

X(R) 

4 

560.0 

-605.0 

81.3 

5.402 

6 

317.0 

-537.0 

82.6 

5.303 

10 

140.8 

-3S3.9 

80.0 

5.31 3 

16 

60.8 

-253.0 

80.2 

5.299 

20 

41.2 

-202.6 

82.2 

5.36  3 

26 

26.26 

-152.8 

79.5 

5.432 

30 

20.00 

-128.74 

80.6 

5.293 

36 

15.42 

-102.46 

80.5 

5.383 

40 

13.21 

-88.23 

79.0 

5.277 

50 

10.33 

-60.10 

78.2 

5.329 

60 

9.56 

-37.91 

77,5 

5.419 

70 

9.68 

-17.17 

77,9 

5.384 

80 

11.40 

6.50 

81.9 

5.31  1 
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Table  A. 2 

INDUCTANCE  AND  CAPACITANCE  OF  THE  CELL 


Vertical 

Horizontal 

Height  of 
Sample 
(inch) 

L3  (MH) 

C4  (PF) 

Height  of 
Sample 
(inch) 

L3  (MH) 

C4  (PF) 

3.564 

155 

8.9 

4.142 

93 

6.4 

3.144 

145 

9.1 

3.292 

93 

6.i.' 

2.374 

145 

9.3 

2.627 

94 

6.4 

1.969 

145 

9.3 

1.9)7 

95 

6.4 

1.229 

145 

9.9 

1.517 

94 

6.4 
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Resistance,  (-J'-) 


LISTING  OF  THE  COKPUTER  PROGRAi-:  WHICH  COMPUTES  THE  L 

CONSTANT  AND  CONDUCTIVITY  OF  THE  SA.VPLE  FROY.  THE  YEA 
lYPEDENCE  AT  ANY  REQUIRED  FREQ'JET.'CY. 

1 0  D IN  F[ 20  ] ’ M[ 20 ] ?  E [ 20  ] • C [ 20  ] ’ DC  20 ] 

20  DIM  E[20]'G[20].H[  20  3.L[S]';:C5] 

20  D I M  PC  20  ] ' UC  20  ]  <  V [ 20  ] . 0[ 20  3 

•^0  DIM  R[  C"  20  3«  ZC  j>  20  3' VC  20  3 

LO  DIM  M[  20  3>Mt203 

t  0  D I M  ijC  20  3  >  JC  20  3  -  U  C  20  3 

70  INPUT  M9-L2 

:E0  IF  M9  =  2  THEN  120 

20  PRINT  "  HOPIZONThL" 

100  PRINT 

110  PRINT 

120  GOTO  lEO 

120  PRINT  "  '"‘ERTIChL  " 

140  PRINT 
ICO  PRINT 
lEO  N5=i9 

170  PRINT  "  L2=".L2 

150  IF  M5=2  THEN  250 

1 20  UC 1 3=2. 325 

200  ::C  2  3=2.475 

210  ::C33=1.31 

220  :'C4  3=1.  1 

2  20  ::l5  3=0.7 

240  GOTO  200 

2‘^0  UC  13  =  2.555 

250  UC2  3=2.5'5 

2-0  UC  2  3=1.305 

2:20  ;:C4  3=1.4 

250  UC5  3  =  0.£5 

200  FOP  1=1  TO  N5 

210  FEtiD  F[  I  3»  hC  I  3 

320  EC  I  3=1 . 3^2-F I-FC  I  3-'  lE-0  2  • 

2  20  CC  I  3  =  0 

240  DC  I  3  =  -l  •  '■  1 . 3*2-P  I-FC  I  3-''  lE-05  ■  ■ 

250  NEUT  I 

250  DhTH  1 1  0 »  2 ?  0 '  2 ?  0 . 00 1 ' 4 « vO .  00 1  ••  3 0 . 002  >  5  <  0 . 00  25 

2  -  0  D  H  T  M  1 0 )  0 . 0  0  5  j  1 5  »  O' .  0  O'  5 "  2 '!'  >  0 .  0  O'  5  <  2  5  >  0 . 0  O'  5  <  2: 0  <  0 .  0  0  5 

250  DhTR  35 , 0 . 005 n  40  <  0 . 00 ?>  50’ 0. 0 1 ’ 50 • 0. 0 1 2  <  70  <  0 . 0 1 4 

350  IhTPI  50'j  0. 01:5)  50  <  O'.  02'  100^  0. 025 

400  FOR  1=1  TO  10 

410  FEhD  5C I  3 

4  20  ne::t  i 

4  TO  DATA  1,4-1. 4. 1.4.1.  255.  1 . 25:5 

440  DhTR  1 . 035< 1 . 075 . 1 . 055* 1 . 055. 1 . 055 

450  FOR  1=1  TO  H5 

450  RERD  EC  I  3.0C  I  3fGC  I  3.HC  I  3 

4-0  ne::t  i 

4:20  DRTR  5  2. 4 ' -2.  22n  0.  0001 . 0. 01 7 
450  DRTR  52. 25. -2. 2. 0, 0005. 0. 0171 
500  DRTR  5  2 .  205  <  -  2. 2  225 »  0. 0005  •  0.  ‘0172 
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?10 

DhTH 

EE.  -1 

.  4  E  «  O' .  O'  '0  '0  4  '  O' .  O'  E  E 

^  IOi 

tHTfl 

C'  E  .  4  .  — 

1 .  ?  4  <  O' .  0  0  0  E-  4  >  0 .  O' : 

T  --  0 

ItiTri 

E  E .  E  4  E ' 

-  1  .  4  441  <  0. 001  <  0. 

T  4  0 

I'tiTH 

EE.  1  - -1 

.  10'  0 . 0  01  1  E  >  '0 . 0  ^  E 

vt'O 

L'HTh 

E 1 .  c'  - 

0 .  OE  '  0 .  00  144. 0 . 0'4 

^  F  0 

IhTFl 
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c* 

I-iTri 

E 1 . 4  r  <  - 

O' .  E  E  .  O' .  O'  '0  E '  O' .  17 

'r  u 

THTR 

E  1 . :  4 .  - 

O' .  E  .  O' .  0  0’  E  E  E  •  0 .  E 

'?  0 

IH^H 

r.  1  .  ^  r  '  ~ 

0’ .  5  4  '  O' .  O'  O'  E  4 .  O' .  E  E  4 

V  0 

L  ri  T  h 

r.  1  ,  4  '  ■■ 

O' .  j .  0 . 0  0  E  4  E .  0  .  E  E' 

H  •  1  * 

ThTH 

EO . E  5 EE 

'  -  0 .  4  c  - 1  4 .  >4 .  O'  O'  E  1  '  '4 

H  i  ij 

0  H  Ol 

E  O' .  E  O'  Z  E 

'  ”  0 .  4  E  E  E  .  O' .  0  O'  4  4  .  'j 

^  4  0 

LhTR 

E  O' .  I  >  E  4  4 

.-0.4  0  4  4  .  O' .  O'  O'  4  E .  O' 

*:  0 

ThTh 

E  0 .  0  O'  E  4 

.  -  O' .  4  E'  4  4 '  O' .  O'  O'  4  4  >  O' 

r-'  r  0 

OhTh 

'54.  E445 

-  O' .  4  4'  ;■  i .  O' .  O'  O'  4 1 .  O' 

7'  fi 

IF 

=1  THEM 

E40' 

4  i  0 

FEITO 

PE  44  0' 

.■-0  FOP  1  =  1  Ti; 
'00  FEhD  LC  I  ] 
'10  he::t  I 


I'Pi'^El  4  4 '  4 . 44 '  ?y. 

44 

■'  10! 

FE4^0FE  EEO 

“  J  " 

TR'E!  1  ??'  14T'  IFIl 

'140"  140 

10’ 

IF  114=1  THEM  440 

'  ■Ll" 

“EETjF'E  O' 40 

■  “*  /1 

FOF  t=l  TO  5 

FOF  1  =  1  TO  MO' 

.*1 

FEE, I'  f:  M.  I  ] 

110! 

•!E:,T  I 

1'.' 

t'E..!  1 

ThTh  1  '  1 '  1 '  1  .  1  !  1 

1  •  1 '  1 '  1 

1  '  1  '  1 .  1  '  1 .  E  .  1 

IRTH  1 '  1  ■  1 '  1  -  1  '  1  ' 

1 '  1 '  1 '  I ' 

1  '  1  '  I .  ! .  1 .  E  .  1 

4  0' 

I'RTm  I'l'l'l'l'I' 

1 '  1  •  1 '  1  . 

I  •  1  •  1 .  2  ■  1 .  E  .  1 

■  0 

TR^R  I'l'l'l'l'I' 

1 '  1 '  1 '  1 ' 

1 '  1 '  1 .  1  ■  1 .  E  .  1 

0 

I'RTR  I'l'l'l'l'I' 

1 '  1 '  1 '  1 ' 

1  '  1 '  1 .  1 '  1  .  E  .  1 

T'O 

'^EE-T'I'FE  1040 

■LlI 

I'  R  T  R  O' .  5  ’  O' .  j '  O' .  0  ' 

O' .  0 '  'I' .  0  • 

0 . 0 '  R .  '  R .  ,  M 

r  u 

DhTR  E.'? 

0  0 

D  R  T  R  0.0' 0 . 0 '  O' .  0  ■ 

0 . 0 ' 0 .  0 ' 

0 ,  0  •  0 .  R .  .  R 

10 

D  R  T  fl  4 . 3 

JLO 

I R  T  R  0.5' 0 . 0 ■ 0 . 0 ' 

O' .  O' 0 . 5 ' 

0 . 0  •  O' .  R  '0.4'  R 

I'RTR  E 

4  0i 

IRTR  0.  0 '  O'.  0 ' 0.  0  ' 

O' .  0  -'  O' .  E ' 

O' .  4 .  O' .  5 '  1  .  1  .  E 

IRTR  :4 

C  0 

T'RTh  0,  O' 0.  O' 0, 0  ' 

0 . 0 ' 0 . 0 ' 

O' .  5  '  R  .  R  '  Oi .  7  .  M 

70 

IRTR  E ' 3 '  1  0 

J  .  ;  . 

4-  :  .  " 
4  .  1  .  f  ' 


1  .  -  .  1 


'■0^0  IF  114=1  Them  lO’O'O' 
■440  FEE  TO  PE  IOjO 
1000  FOR  1  =  1  TO  ‘5 
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1 0  i  0 

PE  ':'  T[  1  2 

1 0  -  0 

NE,;-'  I 

1  0  2-:  0 

D  'i  T  H  4.  ^  4  ,  H , . . 

1  0  4  0 

P'  E  :■  T  i_i  F  E  114  0 

1  ozo 

I'PiTm  :2:  .'?<  ?.  1 

10-' 0 

:f  r',.4=i  THE.N 

1 0^0 

P'  E  :  T  LI  P  E  1  j:  4 

1  0  2-  0 

r  2  0'  }  -  ]  TO! 

1^1  "C  n: 

!  I  •?  0  P  t  D  Z I  ‘  •  I  2 
1  1  :•?  ZCI  -  I  ]  =  :[  i  <1  2- i  £-*-.- 
;i_0  me.:t  I 
1120  HE  :t  f 

1  1  0  I'  h  T  ri  0 ' 0 '  0  •  2> .  5  •  1  .  Z  •  -) .  Z .  ^  .  '■• .  Z  .  1  ! 

llZ'j  DmTh  Z5' 30<  2.0«  2-5.  100' Z4.:  • -JZ‘2' 

1  1  E  0  L  P  T  (i  0 ' 0 ' 0 ' 0 . 5 '  1 .  Z '  -4 .  Z .  Z '  ■? .  Z '  ;  P’  •  1 1 

11"0  DhTR  ZZ  '  I-'O- :  O' 2  Z  .  1  00' Z40  >  4Z0 

1  :■  0  iifi'^0  0 ' 0 ' 0 ' 0 .  z '  1 .  z '  4 '  z .  z '  pTz '  : m >  1 

1  1  r 0  0 M ~ M  Z Z '  2! 0 '  2. 0 '  2' Z '  1 0 O'  Z 4 0 ' 4  ■' H 

1  Z  00  Z 4  T  H  0 ' 0 '  0 '  0 ,  Z '  1 .  Z '  4 '  z .  Z '  ? T Z '  1  '  1 ' 

IZlO  D4Th  ZZ' 2:0' 20' 2Z' 100' Z40' 4ii‘ 

1 Z  ZO  Dm T M  0 ' 0 ' 0 • 0 . Z '  1  . Z «  4 ' Z , Z • zTz ' ; H •  1 1 
1  ...  0  D  M  *  M  Z  Z  '  2:  O' '  2'  '0  '  2-  Z  '  1  0’  O' '  Z  4  >  4  i"' 

1  4  0  Li  4  T  rl  01  I  01  '  O' '  Ol  .  ,  1  C-  .1  c  , 


1  O' '  !  0 .  :  '  1  Z '  : 

1  0 '  1 0 .  '  1  z  '  r 

1 0 '  1 0 .  Z  '  1  z '  1  ■ 


1-4  0 

0 ' 

0 '  O' '  0 .  z 

i:zoi 

DhTh 

;:C 

'  50'.  ;  O' O’ 

2  -  P  0 

DmTm 

C ' 

0 ' 0 ' 0 .  z 

1  ^  7  0 

D'mTm 

*  r 

»  7  C'  •  1 1  •  J 

1  zzo 

DmT  m 

p  * 

0  •  0 ' 0 .  z 

1Z2'0 

Dm'’'Pi 

4  0 

^  T*  0  < :  2  0 

1  2  O'  0 

Dm  ”n 

0 » 

0’ '  C’ '  Ol .  z 

1  ;■  1  0 

I'M  7  m 

4  0 

*  iv*  1  1  0 

1  ]  I'J 

0 ' 

O' '  O' '  Ol .  z 

1  ;  2  0 

I  mTm 

ZO 

1  <  1  f;  M 

1  240 

FOP:  1 

TO  Z 

•  :  >:  0  •  ..  ;  •  4  =  .j 

.  Z '  1 .  Z '  4 '  Z .  •  0 

-O'  .L  Z  O'  >  2  “  Oi ,  s  O'  >0 


V  •  - 1  0  '  t  .,1  ..1 

Z .  Z '  'O .  Z '  ]  O' 

0 '  Z  O'  O' «  .j!  O' 

'r  r .  -  c  2  - 


1  0  .  Z  '  1  Z  '  1  Z 
1  O' .  Z '  1  Z '  1  z 
10. Z' IZ'  IZ 
:o. f .  li. iv 


0 '  u 

■i  c  .  <  cr 


4  '  Z  .  Z  '  0  , 
4  Z  O' '  2  O'.  O' ' 


1 2Z0  l;=::c  i  ]-lz 

lO-'O  L4  =  ':[I-n-LZ 

1  22’0  IF  hE2';.2:'  O.Ol  thei,  i4:.ri 

1  -  '-'I  iF  m£:2'  L4  0.01  Then  I4E0i 

IZ'-^O  LZ  =  LZ-''L4 

1400  IF  LZ  0  then  14-0 

1410  NE.’.'T  I 

1  4  Z  0  G  oi  T  Oi  1 Z  Z  0 

14  2:0  fZ=I 

2440  ^ E=J 

14Z0  GOTO  IZ.2'0 

14E0  f:Z=I-l 

1470  tE=I-l 

1  4  2  0  G  ij  T  0  1  Z Ol 

14'40  iz=i 

IZOO  KE=I-1 


l51‘j  GOTO  150G< 

1020  IF  L2  0  THEM  10 00 
1020  1.0  =  4 
1040  i:f.=o 

100  0  G  Cl  T  ij  1 0  'i  0 
1000  10=1 
lO'O  1.0=2 

1  02  0  FOF'  1  =  1  TO  MO 
10 1^0  IF  1  0  =  .i  0  THEM  1020 

1000  ::o=  lo-ol^o:’  ■  ::[  >  o  ]-:.l  i  o  ] 

101  0  ij T  1 0  H 
1 0  2  0  i !  0  =  0 

1  0  0  F [  I  1  =  F' [  1  O'  I  ]  +  '  H [  t  .f; .  T  ] - p [  f  s .  I  ]  ‘ ^ 
1  0  4  0  IF  110  =  2  THEM  1 0  2  tj 
1  000  F0=:[  TO  ]+  ■  2[  1  0  ]-2[  10]-.  0 
1000  GOTO  1020 

1070  F0  =  2[  1:5  +  0  ]+' 211  O  +  O  3-2I  1  0  +  0  ] '-MO 

1020  Or  =  LC  1:0  ]+’ LI  1  0]-L[»  0]  0 

1020  ui  I  ]=Z[t  O'  I  ]+'  zn  O'  I  ]-Z[  1  O' :  ]  7 

1700  W7  =  T[  1.0  J+'  TL  1  0  ]-T[  I  0  ]  -  -TO 

IT' 10  MEOT  I 

iZZO  M2=0 

I  7 2 O'  I  MF'-.IT  F  V  '  :  !0'  •  To 

1740  IF  F2=0  them  2220 

;~00  N2=M2+1 

1700  FOP  1=1  TO  MO 

1770  IF  FI  I ]  =  F2  them  1720 

1720  ME.iT  I 

17  20  T,2=I 

1  2  Cl  0  7  0  =  E  7  F  G  C 1  2''  I  '  +  E  7  F  ■  -  G  [  I  2  1  ■' 

12  10  7  0  =  E  7  F' '.  ij  C  1  E  .  F' '  -  'j  [  1  2  1  ■ 
l2'20  7 1  =VO-C Om HI  1  '2]  ' 

1  2  20  Vl=:iO-2IM'  HI  1  2;  ]  . 

1240  ::2  =  70-C02  ■  HE  I  2  1  ■ 
iCO'O  V2  =  V0-2IM’  Hi  1  >■'  3  ' 

If'Oo  7  2  =  +  7 2+ V2 -  I  2 

1270  T0='  "1+:  :Z  +  'T1-  I'Z  ■-  7' 2 
12’20  T0=  '  -7 1 -'i'Z  +  7 1 :  2  72 

1  C'  2  0  G  Cl  2'  0  El  2 1 0  0 

1200  1=1:2 

1210  ::i=ui  I  ]-Ft  I  ] 

1  220  72  =  MtI]-F[I] 

12  20  'r'2  =  M[  I  3-'  2>F  IfF[  I  3-07- •  lE-02.' .> 

1  9  4  u  Z  6  ”  i  i  !■'!  w  '1'  j:'  I  Z 
1200  G1=0II] 
rrOO  G2=72  ZO 
1270  G2=G2-G1 
1220  R  2=1  1:3 

1220  Ci2  =  ~'i'2.''  ^■'FI-FI  I  3-ZO-'  lE-00  '  ' 

2000  04  =  C  2-TJ7 
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2010  G'?=b  1 » 1  E  +  Oi* ) 

2020  IF  M9=l  THEH  2070 
2020  Q[H2]  =  FCn 

2040  JC  H3  ]  =  C4ra2  +  0.  SE'?  • -  iFE-O. 

2050  l'J[  N?  ]=(L2  +  0. 5c9  >-■  ■;  Fb-4.  ■?77 

2060  GOTO  2100 

2070  0CN3]=F[n 

2080  JC  N2  ]  =  C4:-2. 7005  FG 

20:-'0  liC  M3  3  =  0 .  326 14.  '  c  F3*F6- 1 E-04 

2100  WRITE  <  15<  21  10  ■  M3- FC  I  3,  J[(J 

2110  FORMRT  2F 1 0 . 0 » 2 F 1 5. 3 

2120  GOTO  1730 

2130  3T0P 

2140  EMD 

2150  Tl=H[t3  3-::8 

2160  T2  =  E:CK3  3-TS 

2170  Dl=CCI.S3-Tl-riC^  33-T2 

2180  D2  =  l:CFS  3-T2  +  DCf::5  3^TI 

2130  T1=;:S-FIC  I  3  3-CC  t.S  3 

2200  T2  =  'i£-E[  f  8  3-ri[  KS  3 

2210  T3=T1^T1+T2-T2 

2220  P3=‘..Iil*Tl+D2-T2  ,'-  T3 

22 30  0.3  =  *  -r  1  -^T2+r2*T  1  ;•  -  T 3 

2240  >;i=EC  ra  3^T5-0[  K  8  3-T6 

2250  VI  =EC  f  3  3-T6  +  OC  K8  3-^75 

2260  :'1=:81-P8 

2270  Vl=Vl-08 

2230  M2=P8-T5-03-T6 

2  2  3  0  '  i'  2  =  F'  8:  -■  T  6  +  C  3  -  T  5 

2300  :82=.:2-E[t  3  3 

2310  V2  =  V2-0CI  8  3 

2320  ::4=::2^::2+v2-v2 

2330  M 3=  < ; 1  *!  !2‘*-'.' 1  *-'i'2  '  74 

2340  V3=':-Kl-V2+Vl-72  •  ::4 

2350  MCt'3  3=ECt  8]*::3-0C(  8  3-V3 

2  360  MC  F3  3  =  EC  f  S  3- .'O+OC  f  3  3^.: 3 

2370  RETURN 

2  380  I MPUT  J 1 »  J2 - J 3 «  j4 - J5 • J6 

2330  INPUT  N4 

2400  IF  N4=2  THEN  2510 

2410  SCfiLE  0.2jJ1.J2 

2420  Mm;! IS  J1>1 

2430  VhMIS'  0'  J3 

2440  FOR  1=1  TO  N3 

2450  F2=L0G':C!C  I  3  ' 

2460  F3=JCI3 
2470  PLOT  F2-F3 
2  4  8'  0  PLOT  F  2 F  3  .<  2 
2430  MENT  I 
2  5  0  0  G  0  T  0  2  3  3  0 
2510  SCALE  0-2JJ4.J5 
2520  :;ANIS  J4,1 
2530  VRNIS  0<J5 
2540  FOR  1=1  TO  N3 
2550  F2=L0G(Q[  I  3.' 

2560  F3=U[ I  3 
2570  PLOT  F2»F3 
2580  PLOT  r2-F3.2 
2530  NENT  I 
2600  GOTO  2330 
2610  END 


4  4  U  4  o  ‘ 

■»F: 3'^  1 E-04  ' 


3  3-  MC  N  3  3 
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r.ibic  1 


nii-r  nRr.NT  rr.Asim.r,  oi'timiziu')  Moi'tn  pauamiih  rs 


No. 

Soil  Type 

a 

b 

c 

d 

c 

r 

s 

-4 

k  X  10 
(mho/cm) 

1  (a) 

Illite 

.‘>(,'4  0 

.00  55 

.3596 

.5974 

59.9 

7‘) 

29 

(b) 

.0058 

.2903 

.4792 

■>8.5 

79 

54 

2  (a) 

Illite 

.5885 

.0060 

.2834 

.5184 

21.8 

79 

16 

(b) 

.5078 

.0027 

.1965 

.6129 

42.9 

79 

26 

(c) 

.4916 

.0052 

.1836 

.6265 

47.4 

79 

28 

3  (a) 

Yolo  Loam 

.5376 

.0053 

.1873 

.6516 

26.1 

79 

21 

(b) 

.5820 

.0061 

.2252 

.6130 

18.5 

79 

17 

Table  2 

COMPARISON  OF  WATER  CONTENTS 
TO  THE  RESULTS  GIVEN  IN 

CORRESPONDING 

TABLE  2.1 

Water 

Water 

No. 

Soil  Type 

Content  (%) 
-Calculated 

Content  (%) 
-Measured 

1  (a) 

Illite 

.904 

2.887 

157.9 

1  59.0 

(b) 

2.592 

3.904 

236.2 

159.0 

2  (a) 

Illite 

0.301 

1.209 

54.9 

60.0 

(b) 

1 .068 

1.506 

86.4 

60.0 

(r) 

1.557 

1.405 

100.4 

60.0 

3  (a) 

Yolo  Loam 

0.408 

0.755 

4  2.2 

4  5.0 

(b) 

0.250 

0.777 

26.6 

45.0 
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Table  3 


OPTlMIZlin  MOHEL  PARAMETERS  TO  PREDICT 
WATER  CONTENT  OF  THE  SOILS 


• 

o 

Z 

Soil  Type 

a 

b 

c 

d 

s 

-ii 

X  10 

(mho/cm) 

1 

lllite 

.0033 

.3396 

.3974 

39.9 

79 

29 

2 

lllite 

.7130 

.0032 

.2376 

.6386 

30.7 

79 

17 

1 

llhto 

.r,2'(9 

.0070 

.2307 

.3988 

2  3.9 

79 

19 

4 

lllite 

.0049 

.2593 

.3929 

1  3.8 

79 

17 

5 

K-lllite 

.  sr.^,  1 

.0047 

.2194 

.6123 

33.8 

79 

21 

f. 

Li-llllte 

.3883 

O 

o 

.2834 

.3184 

21.8 

79 

16 

7 

lllite 

.3614 

.0064 

.3254 

.4203 

22.6 

79 

1  1 

8 

Kaolinite 

.5869 

.0044 

.3327 

.4332 

18.7 

79 

3 

9 

Yolo  Loam 

.5376 

.0053 

.1873 

.6516 

26.1 

79 

21 

10 

Kaolinite 

.6678 

.0001 

.2788 

.5824 

9.8 

79 

78 

11 

Kaolinite 

.8693 

.0261 

.2542 

.7076 

12.3 

79 

6 

T  able 

COMPARISON  OF  WATER  CONTENTS  CORRESPONDING 

TO  THE  RESULTS  GIVEN  IN  TABLE  2.3 

No. 

Soil  Type 

c 

p 

Water 

Content  (%) 
-Calculated 

Water 

Content  (%) 
-Measured 

1 

lllite 

.904 

2.887 

137.8 

159 

2 

lllite 

.542 

.873 

51.5 

54 

3 

lllite 

.353 

.906 

45.8 

42 

4 

lllite 

.179 

.809 

33.9 

32 

5 

K-lllite 

.724 

1 .09  1 

6  3.9 

61 

6 

Li-lllite 

.301 

1.209 

34.9 

60 

7 

lllite 

.322 

1.823 

77.9 

77 

8 

Kaolinite 

.233 

1.616 

67.3 

60 

9 

Yolo  Loam 

.408 

.753 

42.2 

4  3 

10 

Kaolinite 

.772 

.077 

32.2 

29 

11 

Kaolinite 

.461 

.1 16 

21.9 

25 
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Tabic  5 

THE  VALUES  OF  VERTICAL  AND  HORIZONTAL  FORMATION 
FACTORS  AT  DIFFERENT  VOID  RATIOS 


MARYSVILLE 

RED  SOIL 

No. 

e 

n 

•^H 

F 

Log  F 

-LOG  n 

1 

2.1S7 

.686  5 

2.409 

1.850 

1.14 

2.04 

.309 

.164 

2 

1.753 

.6363 

2.695 

2.07 

1.14 

2.28 

.358 

.196 

3 

1.579 

.6124 

2.781 

2.09 

1.15 

2.32 

.366 

.213 

4 

1.41S 

.5868 

3.056 

2.23 

1.17 

2.51 

.399 

.2  5.; 

5 

1.239 

.5534 

3.320 

2.40 

1.17 

2.71 

.432 

.257 
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Table  6  ,  . 

i 

THE  VALUES  OF  VERTICAL  ANO  HORIZONTAL  FORMATION  \ 

FACTORS  AT  OIFFERENT  VOID  RATIOS 


ILLITE  (100%) 


No, 

€ 

n 

^V 

f^H 

II 

F 

Log  F 

-LOG  n 

1 

1.182 

.542 

2.84 

2.67 

1.03 

2.73 

.4  30 

.260 

2 

1,066 

,514 

3.19 

2.82 

1.06 

2.94 

.469 

.289 

3 

.9141* 

,486 

3,42 

2.92 

1.08 

3.09 

.489 

.314 

k 

.821 

.451 

4,00 

3.40 

1.08 

3.60 

.556 

.346 

Table  7 

THE  VALUES  OP  VERTICAL  ANH  HORIZONTAL  FORMATION 
FACTORS  AT  DIFFERENT  VOID  RATIOS 


SNOW  CAL  (95%)  ♦ 

BENTONITE  (5%) 

No. 

e 

n 

'^H 

•^H 

F 

Log  F 

-LOG  n 

1 

1.594 

.614 

2.17 

1.90 

1.07 

1.98 

.299 

.212 

2 

1.458 

.593 

2.34 

2.07 

1.06 

2.16 

.335 

.227 

3 

1.350 

.575 

2.59 

2.17 

1.09 

2.31 

.364 

.241 

4 

1.234 

.552 

2.71 

2.28 

1.09 

2.42 

.385 

.2  58 

5 

1.123 

.529 

2.69 

2.28 

1.09 

2.42 

.383 

.277 

6 

0.998 

.499 

3.1 1 

2.53 

1.10 

2.72 

.435 

.302 

7 

.0778 

.438 

3.37 

3.06 

1.10 

3.16 

.499 

.359 
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Tabic  8 

THE  VALUES  OP  VERTICAL  ANH  HORIZONTAL  FORMATION 
FACTORS  AT  Oil  PERENT  VOlO  RATIOS 


SNOW  CAL  (60%) 

♦  ILLITE  (40%) 

a 

No. 

e 

n 

^H 

1-^ 

II 

< 

F 

Log  P 

-LOG  n 

1 

1.4^0 

.589 

2.68 

2.29 

1. 08 

2.42 

.384 

.230 

2 

.494 

3.41 

2.82 

1. 09 

3.02 

.479 

.302 

.8^4 

.455 

3.59 

3.25 

1.05 

3.36 

.527 

.342 

.748 

.428 

4.34 

3.72 

1.08 

3.93 

.594 

.369 
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OPTIMIZED  MODEL  PARAMETERS  AT  DIFFERENT  WATER  CONTENTS  -  SNOW  CAL  (60%)  +  ILLITE  (40%) 


Cr  o  ro  ov  CO  rs4 

^  tT'  kr\  — ■  — . 

tr\  ct  Csi 


cj-  OO  vO  ON 

^s.  vr\ 

O  ^  ^  vD 


s£>  <y\  vO 

tTN  \C  v£>  oo  O 


OO  On  W-N  sC 

—  r>.  \£>  ir\  \0 
oo  vXi  -3* 


pv.  V.D  sD  (Nj 
V.0  vO  vD  vO 


^  O  OO 
<N  O  00 


•3“  ^©0^<^  — 
cP  r*^  fsj 


00  r*>.  <JN 

^  OO  CM  O  00  >3 

—  O  CTv 

0^  oO  oo 


\0  0^  rv.  vC  ^  oo 

s£?  vD  vD 

^  v0  pN.  fN^ 

o  o  o  o  o  o 


m  CM  po  pv  r^J  K. 

Px  (Js  ^  cj- 

o  o  o  o  o  o 


fN  sC  ^  —  CsJ  (T' 
vO  —  O 

— .  O  ^  J-  \D 

ps.  p>.  p*,.  p^  Pv  px. 


-rp  oo  —  'X  ct  O 

CN  -3  P^  oO  r^J 

—  O^  ^P^  —  ^ 

ir\  pn  p*^  fNJ  fNJ 


— •  ri  p^  ^  \r 


28.35  .69U6  .032  .1279  .8157  26.8  16.4  6.5  .426  .335  .749  1.33 

24.54  .f:,478  .082  .1315  .7969  22.3  15.4  6.2  .314  .343  .648  0.94 


OPTIMIZED  MODEL  PARAMETERS  AT  DIFFERENT  WATER 
CONTENTS  -  SNOW  CAL  (95%)  +  BENTONITE  (5%) 


— •  OO  rj  -O  f  C-O  CT' 
r  I  oO  r  r  i  sD  -  •  '»s 

^  rfy  c<\  Cl 


0<^rvi— .^OOCOOO 


r>.c}-  — ’  ^  ^  ^  'Xi  oo  \0 
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fTx  V  *  ♦ 
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rsl  — •  'T'  oo  oo  vX'  ct 


r'v  ct  ro,  — .  — <  <Ni  ON 
rj  — «  ON  t5  — •  CO  t •>  r>. 


irNiAuN,r^<^rN(S4Csj— » 


OO^Ai^OnO^NOnw^U^ 

^oc^a^oA^^o^v  —  cpo 

t'v'Ar'N  —  ONfNK^CvjfSv 


—  -^oOCvJvDf'^OcJ-ON 
<TvNX;‘ArN.c'NsCrsicv4sc 

hsrNrNC5-w*\tr\<j'<3-  — 

oocooor>oooooooooc 


0<Ar^fM4r\oO<T'r'J^ 
rNO'sO'Aooo^t^-ON— • 
>OOOOOON(7N<7^00— ' 
rNOOOOO—  — *— 


ooo  —  OnCnOt^oO— - 

iTNU^CsIoOcf'A^OOcJ- 

ooooooooo 


fNoocjN  —  — 

for«NsCr^»^  — ' 
rNs-ct  ONfMr^^oocTN  — 
vDvCsCnOnC'^vCnOv£> 


—  ^fv,  —  r^sCr^rxi^ 
sC*~'<T^oO(jNf^f'-'X‘N£) 


-;f  OO  r»''  ON  '/^  — •  rN  ON 

vDVT',  4rN^ct-:^^r^r»^rNl 


— ‘CNif^^AivDr^ooON 


h-  r^r*NrJv^'T' 
O'  I  -  '  -  O  f  r 


pN  fo  r<N  (Nj  r  J  r -J 


NC-:s-r>.oo«:orsirv  vr-ic' 
^  O''  ^  r  i  o  CO  ^'- 

r><  lA  ^  04  —  VO  '  - 


IT,  rc  ir\  ^  r  J  cj  >-'  ..-« 

r<'  r«N,  f»A  r'N  r*"',  '••  .  .  -v" 


-  I  I  -IvX'-^)'- 
ON  lA  CN  ln  O 
rA  r  4  — ^  (.'N  ON  CO 


ON  On  — -  X 
I  -  t  .  -  X  X'  X  X 


OvONON^tONiT'  —  rsirsi 
'A  —  oor'N.iAf^  C*  — 


COO^O^Ac^onO'  — 
COONiA^ooNsTO  'y- 

rv^u-vf*N—  ONh^::?.  OnC 

--a; 


ONO'^i^rNvrN'vr^N  ^ 

—  O  nX  rA  po  00  O  Aj 
tNs.hv.UA»/N»A4S-0'  CN 
COOOOOoOoOOOSOCO 


NTTNlONTrsiAl—  C^F's 
—  —  ON^}-^  r-.vX 

&OCOONO^o^o^c7^  O  — 

oooooo— - 


ctrA^ONirsiAsXNXA,' 
cONTr^cj-  —  c^— •rN.c*' 
OOOOOOO  r'O 


rsi  —  rNjo''pNi/~.,^  u-'iT' 

rNv>ocf-«:?-r'.rsjr<'  c^  — 
pN^iArNipNCM  —  C'O- 
XnC'sTnXvCnCX  XA 


—  sx  —  ^  A  o 

—  ro  fvj  — «  p^  j/N  :»; 


vX  O  ‘A  —  vX  Ag  r- . 

X  nC  A  »rN  45  '• 


—  AjrA^iA\Xr% 
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OPTIMIZED  MODEL  PARAMETERS  AT  DIFFERENT  WATER  CONTENTS  -  YOLO  LOAM 


26.61  .7291  .047  .2770  .6201  10.36  10.34  6.32  .036  .663  .731  .129 
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Table 

THE  VALUES  OF  'M*  AND  THE  ELEC'TRK'AL 
INHEX  FOR  nil  I  EUI.NT  SOILS 


Soil  Type 

M 

(Obtained  for 
Cam-Clay) 

M 

(Obtained  for 
l3oiinding  Surface) 

A 

f 

A^/f 

Illite 

2.55 

1.S5 

1.07 

1.60 

.716 

Snow  Cal  . 

Elite  (40%) 

2.32 

l.SO 

l.DS 

1.6Q 

.729 

Marvsville  Red 

Soil 

1.S5 

1.(32 

1.15 

1.73 

.'63 

Snow  Cal  (95%)  + 

P.r»f  ^  t  rsn  i  to 

1.19 

1.45 

1.10 

1.43 

.547 
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Table  14 

THE  UNDRAINED  SHEAR  STRENGTH  OF  NORMALLY  C:ONSOLII'IATEn  SOILS 
PREniCTED  BY  CAM-CLAY  THEORY  AND  BOUNDING  SURFACE  THEORY 


Soil  Type 

q  /p 

'  o 

By  Cam-Clay 

Theory 

By  Bounding 
Surface  Theory 

lllitc 

1.11 

O.Sl 

500%!.  Cal  {9')'\0  + 
Bentonite  (yV) 

0.49 

0.3S 

Snow  Cal  (60%)  + 
lllitc  (40%) 

1.03 

0.90 

Marysville  Red 

Soil 

0.76 

0.93 
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FIG.  2  Variation  of  Dielectric  Constant  and  Conductivity  as  a  Function  of  Frequency  for  Cohesive  Soils 


1  2  3 


O  Solid  Particles 
1^/1  Solution 


Three  components  of  the  current  paths 
through  the  soil  sample. 
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[Z3  Solution 
i — I  Solid 


Figure  Three  Element  Electrical  Model 

Each  Zone  is  Represented  by  a  Parallel  Circuit  of 
Resistor  and  Capacitor  .  Impedance  of  Each  Zone  is 
Determined  by  its  Dimension  (a,b,c,d)  and  Specific 
Conductivity  and  Dielectric  Constant  of  Material  Forming 
the  Zone  (k. .  k, .  ,  «,) . 
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termi:iation  of  electrical  proDcrtien  at  different  ',va'*‘er  eo>n  +  ov^ 


Fig.  9  T;ie  imreiar.ct:  ana].yzi..'r  uaod  Tor  thii 
rr.easurement  oi  elt'Ctricai  properties 
in  th.e  radio  Irequoncy  rcuige. 
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different  water  contents 
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'iD-Ai31  6bO  IN  SITU  CHARACTERIZATION  OF  SOILS  FOR  PREDICTION  OF  J/t 

STRESS  STRAIN  RELATIO..(U)  CALIFORNIA  UNIV  DAVIS  DEPT 
OF  CIVIL  ENGINEERING  K  ARULANANDAN  ET  AL.  10  NOV  82 
UNCLASSIFIED  AFOSR - TR - 83 - 0680  Af OSR- 8 1  - 02 16  F/G  8/13 
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Dielectric  Constant, 
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Dielectric  Constajit, 


Conductivity 


Conductivi  ty.lT-X  10  (mho/cm) 


18 


YOL 


i 


Dielectric  Constant 


Frequency (MHZ) 

Fig.  2i  Dielectric  constant  as  a  function  of  frequency 
at  different  water  content. 
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1  10 

Frequency ( MHZ ) 


Fig.  22  Conductivity  as  a  function  of  frequency  at 
different  water  content. 


Conductivity ,  (TXIO”  (mho/cm) 


100 


10 

Frequency (MHZ) 

lelectric  constant  as  a  function  of  frequency 
t  different  water  contents • 


Dielectric  Constant 


Conductivity, (T;*  mho/cm) 


tl 


/' 


Yifr,  31  Tti^  fall!  ti  "r:  nhif  bct,v/eon  t!iO  ratio  cl 
intra  a:;:!  i  .’‘or  clurtcr  void  ratio,  e-j./ 
and  the  tot.d  void  ratio,  e. 
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Ratio  of 

Infra  ond  Inter  Cluster  Void  Ratios,  ei/ep 


Totol  Void  Ratio,  e 


Fig.  32  The  relationship  between  the  ratios  of  intra  a;ii 

inter  cluster  void  ratios  and  the  total  void  ratios. 
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p  *=  Pre-cosolidation  pressure 

0?R=  p„Vp<, 


34  Schematic  Illustration  of  the  definition  of 
X  iK  tM  »  e  and  OCR 


Slope  of  Virgin  Consolidation  Line.X 


Fig.  36  The  relationship  between  the  pressure  and  the 

void  ratio  during  an  isotropic  consolidation  test 
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Log  (X) 


Log  (%ef>)e  =  I.O 


-1.0  -0.5  0  0.5  1.0 


Fig.  37  The  relationship  between  the  slope  of  the  isotropic 
swelling  line,)v  i  and  the  ratios  of  intra  and  inter 
cluster  void  ratios  at  e=1.0 
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FRICTION  IN  SOIL 


c.  Friction  betv;een  roug 
clay  surfaces 


d.  Friction  between  smo- 
clay  surfaces 


Fig.  38  COMPONENTS  OF  SHEAR  RESI STAiNCE(LAMB=  I960 ) 
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Fig,  The  relationship  between  X  loSgP 

q/p 
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Fig.  '*1  The  relationship  between  the  slope  of  the  critical  state  line,M 


|^2  The  relationship  between  q,  p  and  v  for  undrained 
conditions. 
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Bounding  Surface  Model 


SNOW  CAL{60^) +ILLITE(40?5)— p  =  69.3  psi 

-  Predicted  by  bounding  surface  theory 

-  Predicted  by  Cam-Clay  theory 

O  Experimental  results 


omparision  of  experimental  auid  theoritical  stress-strain  ralationship 


Snow  Cal(95%}-t-Bentinite(5%),  Pq  -56.9  psi 


Compari  si  on  of  6xperiinen"tal  and  theoritical  stress-strain  relationships 


Fig,  Comparision  of  experimental  and  theoritical  stress-otrain  relationships. 


lUite,  Pn=56.5  psi 
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